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Abstract 

In order to develop the methods of thermodynamic analysis for the living cell, two 
models of protoplasm microstructure of the living cell in resting state were suggested. 
Both models are based on the assumption that the Ling's cell as a statistical me- 
chanics system is non-ergodic. In the first, Van der Waals model, the protein-protein 
interactions, which form the physical basis for the cell functioning, are considered as 
a interactions of key importance. It is postulated that protein molecules are situated 
in points of some space lattice (the Ling model of a cell) they assemble to aggregates 
at equilibrium state, corresponding to the dead protoplasm. In the second model we 
consider protein conformation at the resting state and conformation changes while 
the cell is passing from the resting state to the equilibrium state (dead protoplasm). 
The investigation of the models and comparison of their characteristics showed that 
the convenient tool to define the energy minimum of the system under consideration 
is a Hamiltonian describing the superfluid Bose gas on protein configuration space. 
Our approach allows us to define the thermodynamic features of the living (at resting 
state) and dead protoplasm in a new way: in the first case the system is characterized 
by the unfolded state of proteins, in the second case proteins are folded and aggre- 
gated. Obtained results prove the applicability of our approaches for thermodynamic 
characteristics of the Ling model of a cell. 
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1 Introduction 



R. Feynman wrote: "All things are made of atoms, and that everything that living things 
do can be understood in terms of the jiggling and wiggling of atoms." (Feynman, 1963). 
To move beyond this assertion, it is necessary to adopt common principles of organization 
of atoms and molecules in living systems. These principles, compatible with the existing 
analytical apparatus of thermodynamics and statistical physics, have been formulated by 
Gilbert Ling (Ling, 2006). The living cell model created by him was used as a starting point 
of his study. 

According to Ling, fundamental properties of the living cell are explained by the single 
physical factor — sorption properties of its proteins. An unfolded (linear) protein molecule 
binding water (multilayer adsorption) and K + (in the presence of Na + ) under the control of 
ATP represents the smallest part (unit) of a living protoplasm which still keeps the main 
physical characteristics of the living cell. Later Matveev (2005) offered to call the unit as a 
physiological atom or physioatom. 

The main physical state of the phisioatom, and accordingly, the cells comprising them, 
is, according to Ling, a resting state. The physical nature of this state determines, on the 
Ling's theory, all forms of biological activity of the cell and therefore the analysis of this 
state is a key issue of the physical theory of the living cell. 

Compatibility of the Ling's resting cell organization principles with analytical meth- 
ods of modern theoretical physics was first shown in our previous work (Prokhorenko and 
Matveev, 2011). The base of our approach is the fact that the majority of statistical me- 
chanics systems (including the most realistic ones) are non-ergodic that was proved by one 
of us (Prokhorenko, 2009). The generalized thermodynamic analysis (generalized thermody- 
namics) of Ling's cell we proposed (Prokhorenko and Matveev, 2011) allowed us to explain 
(in framework of adopted boundary conditions) a number of physiological phenomena that 
occur when cell is in activated (excited) state: exothermicity of transition to excited state, 
change of cell volume, folding of natively unfolded proteins (which determine, by Ling, a 
main features of the resting state), efflux of cell K + and wider — major redistribution of 
physiologically important ions between the cell and its environment. However, we deter- 
mined the sign (direction) only of these processes, there were no numerical evaluations. In 
other words, the results were obtained as inequalities. These are the normal features of 
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thermodynamics (not of the generalized one only): it allows us to obtain general relations 
between thermodynamic variables which are independent on a nature of intermolecular in- 
teractions. It is needed to construct models with desired properties (including microscopic 
models) to obtain specific numerical values of thermodynamic variables and then investigate 
them by methods of statistical mechanics. 

Indeed, inequalities we have determined (Prokhorenko and Matveev, 2011) were based 
on the postulate of relative entropy maximum for the resting cell standing in the state of 
thermodynamic equilibrium with environment. However, an equilibrium with environment 
does not mean the equilibrium state (the absolute maximum of entropy) of a system; that's 
why the relative nature of the entropy maximum is indicated. We consider the resting cell as 
a system in a steady non-equilibrium state described by the generalized Gibbs distribution 
(Prokhorenko and Matveev, 2011). In other words, the case is the maximum among of all 
the states described by generalized Gibbs distributions and constructed using the fixed set of 
first integrals in the involution. These inequalities give evidences of the negative determinacy 
of the second entropy derivative matrix with respect to parameters describing a system in a 
state corresponding to the relative entropy maximum. 

Involvement of statistical mechanics methods brings up an issue of certain properties of 
the investigated system. Ling's model of the living cell gives, in our opinion, an interesting 
material for such analysis. So, after construction of the generalized thermodynamics of 
the resting state, set out in (Prokhorenko and Matveev, 2011), it makes sense to turn to 
some of the structural characteristics of the investigated system (as it often happened in 
the history of thermodynamics and statistical mechanics). In our case, the problem arises 
of constructing various models of protoplasm, and their investigation by various (mostly 
approximate) methods of theoretical physics. In this paper, the authors make the first steps 
in this direction. 

The first model, we call the Van der Waals model, focuses on the nature of interactions 
between protein molecules only. As Ling, we assume that protein molecules in the resting 
cell embedded at specific sites in the lattice of a crystal, and the distances between proteins 
are so long that interaction between them can be neglected. This assumption makes it 
possible to use the method for calculating thermodynamic potentials of ideal systems in 
order to determine thermodynamic potentials of the resting protoplasm with a specified 
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structure. In the case of dead protoplasm (the state opposite to the living resting state), the 
protein molecules are associated due to secondary (non-covalent) bonds in large equilibrium 
aggregates. In this case, to calculate thermodynamic potentials we use the formula, we have 
obtained in this paper, for corrections of free energy values in the case of formation of large 
aggregates. 

As part of the Van der Waals model we have obtained (i) numerical estimation for heat 
amount released when erythrocyte die, (ii) the estimation for number of protein aggregates 
appeared in dead protoplasm, and (iii) the estimation for fraction of whole cell volume 
occupied by these aggregates. All these estimates are in good (qualitative) agreement with 
available experimental data. 

In our second model, we also assume that protein molecules in the resting cell embedded 
at specific sites in the lattice of a crystal, but at this time the focus is on internal structure 
of a protein. We consider this model at zero temperature (the energy scale), which makes 
it possible to use (to calculate the ground state) some effective Hamiltonian describing a 
superfluid Bose gas in the configuration space of a protein molecule. Based on the repre- 
sentations set forth in (Prokhorenko and Matveev, 2011), we define parameters of effective 
Hamiltonian corresponding to living and dead states of protoplasm; we show (in accordance 
with our assumptions) that proteins in the resting state (that determine key properties of 
the system) are natively unfolded, whereas in dead protoplasm the same protein molecules 
are folded. 

In Appendix 1 the process when the living protoplasm transforms into a dead one is con- 
sidered as it appears in the physical point of view. In Appendix 2 we discuss the mechanism 
by which ATP is able to effectively influence sorption properties of proteins of Ling's model 
for water and physiologically important cations. 

2 Non-Ergodicity and Crystallization 

Let's consider some non-trivial issues that arise when we consider the cell as a nonergodic 
system. According to the Ling's model (Ling, 2001), water in the resting cell is in a bound 
quasi- crystalline state (important, water content is about 44 mole/kg wet weight of the cell). 
The bound state of water and its massive amount in the cell has fundamental importance for 
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the understanding of physiological processes (Ling, 1997). Therefore, one of key issues of cell 
physics is the question: which properties the crystal has in terms of our recently proposed 
approach, generalized thermodynamics (Prokhorenko and Matveev, 2011) Let's begin the 
consideration of this issue with finding out the relation between non-ergodicity of a system 
(for example, the Ling's cell) and its solidifying capability at low temperatures. 

Despite this problem definition is non-physiological, its solution will allow us to verify 
once more that a large number of systems of statistical mechanics, including our model, have 
a non-ergodicity property. Our argument will be largely heuristic rather than rigorously 
mathematical character. In the mathematical physics the rigorous theory is often preceded 
by formal theories handling objects of poorly ascertained mathematical context. However, 
the heuristics presented here are of interest, in our view, as a basis for more rigorous methods. 

At first, let's define the ergodicity for statistical mechanics systems. 

Definition. Suppose the quantum system is described by Hamiltonian H and K 1 , K\ 
are some commuting (among themselves) self-adjoint integrals of motion. The system is 
called ergodic with respect to the set of integrals Ki,...,Ki if any dynamical variable com- 
muting with H, Ki, K\ is their function. 

To give a classical analog of this definition we should just replace the word "commutator" 
by a Poisson bracket everywhere. 

Usually, the operator of system momenta P and operator of particles number N are used 
as trivial integrals. 

At first let's show how the non-ergodicity of a system arises from its solidifying capability 
at low temperatures. Let's consider the system at solidifying temperatures, supposing them 
T < T T > 0. In addition, the system can move through the space as a solid body and its 
coordinates (as as solid body) are six real numbers 

X U X 2 , X 3 if i, f 2 , <P3, (!) 

where x±, x 2 , x 3 are Cartesian coordinates of the system's center of mass, and ipi, <p 2 , <^3 
are some coordinates characterizing the position of a system (as a solid body) relative to 
its center of mass, for example, Euler angles. Let p±, ..p 3 ,7Ti, ...,tts be momenta canonically 
conjugated to them. The Hamiltonian of the whole system 

#(xi,...,X3,</?i,...,</?3,Pi,...,P3,7ri,...,7r 3 ) (2) 



5 



is a function of variables xi, x 2 , x 3 ipi, <^ 2 , <p 3 , conjugated momenta to them and operator 
"in other variables". Free energy of a system is given by: 

F(xi,...,X3,^i,...,^3,pi,...,p3,7ri,...,7r 3 |T) 

H(x 1 ,...,x 3 ,<p 1 ,..., V3 ,p 1 ,..., P3 ,n 1 ,...,Tr 3 ) 

= — Tlntre T , (3) 

where trace is taken by Hilbert space which is "left" after separation of variables de- 
scribing the motion of a system as a solid body. We won't refine the meaning of 
words enclosed in quotation marks considering them as intuitive clear. It's clear that 
F(xi, £3, <fi, f 3 ,pi, ...,P3, 7Ti, n 3 \T) does not depend on variable f 3 (if, for example, 
ipi, if 2, <p 3 are Euler angles). 
But 



H(x 1 ,...,X 3 ,ip 1 ,...,ip 3 , Pl ,..., P:j ,TT 1 ,...,TT 3 ) ^— v X i (x 1 ,..,TT 3 ) 

lntre t = ^ ^(xi, 7r 3 )e t , (4) 

i=0 

where di = dimLj is a dimension of eigenspace Li of operator H(xi, 7r 3 ) corresponding to 

00 

Xi eigenvalue of operator H(x± : 7r 3 ). But since ^ o?i(a;i, 7r 3 )e~ Ai ^ 1 '--' 7r3 ^ (/3 := ^) does 

i=0 

not depend on ip 3 (for different then di(xi, ...,7r 3 ) and Aj(xi, ..,^3) are also independent 

00 

of (p 3 . Indeed ^ di(x±, 7r 3 )e _A ^ :El '"' 7r3 ^ is really a Laplace transform of a measure 

00 

y^£(A - Ai(xi, .., 7r 3 ))tZi(a:i, vr 3 ). (5) 

i=0 

So, for all values </? 3 when others values parameters xi,...,tt 3 are the same, operators 
H(xi, x 3 , ip 1 , <ps,pi, ...,P3, 7Ti, 7r 3 ) are unitary equivalent. And after making the 
appropriate unitary transformation of H depending on Xi,...,ir 3 , we can conclude that 
H(xi, x 3 , <fi, <f 3 ,Pi, ...,P3, 7Ti, 7r 3 ) do not depend on </? 3 . We have considered vari- 
ables xi, X3, ipi, tp 3 ,Pi, -.p 3 , 7Ti, 7r 3 as classical ones since they describe macroscopic 
degrees of freedom and appear to be very large. Now we again consider ip 3 , n 3 as quantum 
variables, replacing them by corresponding operators ip 3 , if 3 , and considering the Hamil- 
tonian H i (x 1 , x 3 , ipi, ¥2, Pi, ■■-,P3, ^1, ^2) obtained by replacing the variables ip 3 , n 3 with 
the corresponding quantum- mechanical operators (fi 3 , n 3 , in H to describe our system. This 
operator acts in Hilbert space H ® V where T is a Hilbert space corresponding to operators 
y5 3 , if 3 . The fact of independence of H(x±, x 3 , <f±, f 3 ,Pi, ...,P3, 7Ti, n 3 ) of <p 3 is stated 
now as commutativity of H 1 with tt 3 , and the presence of nontrivial first integral of a system 
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means of course the Hamiltonian degeneracy. The last fact can indicate the non-ergodicity 
of the system, but this new integral should be commutative with the momenta operator. It 
can be achieved for example by consideration of the integral U := ^ instead of 7f 3 , where G 
behaves like ~ V 3 if the volume V of the system approaches infinity. Then II asymptotically 
commutates with the momenta operator which means the non-ergodicity of the system. The 
reason for choosing G ~ V% is shown below. 

However, instead of considering II as a new independent integral we prefer other way. 
Variables x±, x 3 , ip±, </?3,Pi, •••,P3, 7Ti, 7r 3 are canonically conjugated variables satisfy- 
ing to the Hamilton's evolution at temperatures T < T for the Hamiltonian F(xi, vr 3 |T) 
(see section 3). Just as we did before, we can show that F(x±, ...,ir 3 \T) does not depend on 
Xi, £3, (fii, </? 3 . But the system in the equilibrium state has the energy (in the thermody- 
namic limit) proportional to the volume of this system, therefore E < CV for some constant 
C. On the other hand, if ii, I 2 , h are the eigenvalues of the inertia operator of our system as 
a solid body and Ui,u 2 , uj 3 are components of angular velocity along corresponding principal 
axes of inertia operator, then E > + + ywf . However, 7i, I 2) -^3 ~ V%. Therefore, 
in the thermodynamic limit uj\ = u 2 = ^3 = and our system can moves only parallel to 
itself in other words, (p\ = (p 2 = 03 = 0- So, variables ipi, ip 2 , y?3 are motion integrals of the 
system commutating with the impulse operator which makes our system non-ergodic one. 

Since CV >E> f tuf + f u% + f u\ = + + and I U I 2 ,I 3 ~ Vl then 7^. vr 2 , tt 3 ~ 

vl. 

At low angular velocities uji, ...,u; 3 the free energy of the whole system is presented as 

F(x u vr 3 ) = F ( Xl , ...,p 3 ) + M + ... + (6) 

for some function F (xi, ■■■ 1 pz) of variables Xi,...,p 3 . Momentum variables 7Ti, ...,7r 3 can be 
chosen so that time rates of change of canonically conjugated coordinates may be equal to 
uji, ...,003. Then ix\ = iiWi,...,7r 3 = I3U3 and the effective Hamiltonian of the system equals 
to 

2 2 

F( Xl , vr 3 ) = F {x u ...,ps) + + ... + (7) 

The form of this Hamiltonian together with the fact that in the thermodynamic limit 

4 

7Ti- 7T2, 7T3 ~ Vs imply that ipi, ip 2 , V?3 are the integrals of motion commutating with an 
operator of the total system momenta. 
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Now let us ask the question arising out of the context of our analysis: why the crys- 
talline state of matter is stable at temperatures different from zero. The fact that matter 
solidifies at zero temperature is almost evident: configuration of the system must achieve 
the minimum potential energy. The question arises: why do atoms keep on doing just small 
oscillations around points of the lattice at non-zero temperature and why does the lattice re- 
main faultless, though thermal fluctuations seem to break it. This question is closely related 
to the question of why Ling's cell is stable at relatively high temperatures. We shall try to 
answer with help of our generalized thermodynamics (Prokhorenko and Matveev, 2011). 

So, let H be a Hilbert space of our system, H is a Hamiltonian of our system, and E 
is the lowest number belonging to spectrum, and E is a spectral projection of H onto the 
eigensubspace H corresponding to the eigenvalue E Q . 

In classical terms, at T = atoms composing the system have an arrangement which 
meets the condition of minimum potential energy. This means the body solidifies at T = 0. 
As we suppose, at that moment atoms are situated in points of a crystal lattice. But 
according to the accepted approach this lattices are not invariant under infinitesimal rotation, 
i.e. the system state obtained from the initial one by an infinitesimal rotation does not align 
with the initial one. This results degeneracy of E , i.e. trE > 1. 

Now let's complete {H} to obtain the complete set of (commuting) observed values by 
self-adjoint operators Ki,K 2 , .... Here we use Dirac terminology (Dirac, 1958). 

Completeness of the system of observables H, K±, K 2 , ■■■ means that their joint spectrum 
is simple (non-degenerated). Let Pi,P2, ••• be orthogonal projectors in "H projecting to the 
subspaces of subspace lmE (i.e. PiE = E Pi = Pi) and are projections to their own 
subspaces of operators family H, K±, K 2 , .... All P 1 , P 2 , ... are clearly one-dimensional due to 
completeness of operators family H, K ly K 2 , ■■■■ The generalized microcanonical distribution 
(Prokhorenko and Matveev, 2011) describing our system, can be taken, for example, in the 
following form: 

P = Pf, (8) 

for any /. 

Observed values of the integrals Ki, i = 1, 2, ... in the state p are given by the following 
formula: 

K' t =tr(pK t ). (9) 
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Let SO (3) be a group of self-rotations in Euclidean three-dimensional space, o is an arbitrary 
element of this group, 6 is a unitary representation of o in the state space H of our system. 
When subjected to transformation o G SO (3) the state p comes to 6pd + ). Let's show that, 
if necessary, replacing the complete set H, K\, K2, ••• of observables by another complete set 
H, Li, L 2 , ... allows us to choose / from (jSJ so that for some o G 50(3) and for some integer 
i 

tr(pKi) ^ triopo+Ki). (10) 
If the stated conclusion is false, then Vz, j = 1, 2..., Vo G 50(3) , we have 

trioPid+Pj) = triPiPj) = Sij. (11) 
But the latest means that Vi = 1, 2, ... 

oP,d + = P^ (12) 

Let fi be a unitary vector stretching the image Pj. It follows from ( f!2|) that Vz = 1,2... 

oft = exp(i(pi(o))fi (13) 

for some functions (fi(o) on 50(3). The last conclusion is true for any complete set of ob- 
servables H, K\, K2, ■■■■ In particular it is clear that in ( !T3|) we can choose arbitrary the 
orthonormal basis {fi} in lmE . Thus Vo G 50(3), the restriction of 6 on lmE should be 
diagonal in any orthonormal basis, therefore, the restriction of on lmE must be propor- 
tional to identical operator. But 50(3) has no one-dimensional representations except the 
trivial one. Therefore, Vo G 50(3) 6 = 1. Thus, any ground state of our system under each 
rotation should come to itself; but it is false as we seen above. The statement is established. 

So, the generalized microcanonical distribution Pf have the property that some rotation 
of a system causes the change of integral Ki averaged over this state for some f,i = l, 2, .... 

Now, if we give a sufficiently small non-zero temperature to our system, then (as it follows 
from the principle of physical continuity) for this temperature there exists a (generalized) 
equilibrium state described by a generalized microcanonical distribution p, such that after 
some system rotation the corresponding observable value of integral Ki must change for some 
i = 1,2.... But the system entropy does not change under the rotation of the system. This 
means that for a fixed energy (corresponding to enough small temperatures) the system 
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entropy has a plateau of d > dimensions that provides the stability of our generalized 
microcanonical distributions, as it was discussed in our previous work (Prokhorenko and 
Matveev, 2011). 

In addition, let's remark that instead of speaking about the completeness of the sys- 
tem of observables H, Ki, K 2 , we should speak about the macroscopical completeness 
of this system. We say that the system of macroscopical observables quantities Oi, 2 , ... 
is macroscopically complete if any macroscopical quantity (which clearly commutates with 
Oi, 2 , ... because all macroscopical quantities are simultaneously measurable) is a function 
of observables quantities d±, (9 2 , .... 

Thus, for enough small temperatures the system have stable stationary states which are 
described by generalized microcanonical distributions which are not reduced to the common 
microcanonical distribution. We identify the crystal states of the matter exactly with such 
states, and the stability of the crystal state could be explained by just proven stability of 
the corresponding microcanonical distribution. 

3 Van der Waals Model of Protoplasm 

The general description of this model is presented in the introduction. Within the framework 
of this model and basing on our generalized thermodynamics, we give numerical estimates 
for some changes proceeding in a cell while it excites or becomes damaged: the amount of 
emitted heat by the cell and amount of released potassium ions from the cell to environment 
(according to Ling, potassium ions in the resting state are bounded by proteins). 

Let's consider two extreme protoplasm states: resting state and "dead" protoplasm. 
First, let's discuss how does the dead protoplasm appear in the context of our model. 

In the dead protoplasm protein molecules are in the folded state Prokhorenko and 
Matveev, 2011), and we suppose they are homogeneous balls of radius r and dielectric 
permittivity e'. The dielectric permittivity of the other matter in the cell is denoted by e. 
Assume that M is a mass of a protein molecule. Since a protein molecule contains a lot of 
atoms, M is a very large quantity. 

Let's write out the equation of motion which describes the motion of protein molecules. 
We suppose that the cell is described by classical mechanics, however further obtained results 
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clearly imply that the answer for the quantum case is the same. The motion of protein 
molecules may be considered as classical motion because of a high value of M. Let's denote 
x = (p, q) are coordinates and momenta of all protein molecules. Put by definition that 
V = ip'iQ 1 ) are coordinates and momenta of all other protoplasm components. Let H(x,y) 
be the Hamiltonian of the whole protoplasm. Then Hamilton's differential equations on x 
take the following form: 

. = dH(x,y) 
dq 

i=«agJ!>. (14, 

However, since the cell is dead, its distribution function is a Gibbs distribution function. 
In particular, the conditional probability density that variable x takes value x' provided 
variable y takes value y' is given by: 

w(x \y) = Zi( y) e T > ( 15 ) 
where T is a system temperature and 

Z x (y') := / dxe-^r 12 . (16) 



Since protein molecules move slowly and their mass M is very large (thousands of D), in 
(|14p we can replace right parts by their averages over distribution w(y\x). Omitting rather 
trivial calculations, we find that the averaged system ( Il4p is a Hamiltonian one too and the 
corresponding Hamiltonian is a free energy of the system. More specifically, the averaged 
system (fll]) is given by: 

dF(x\T) 
dq 

« = ™ < 17 > 



V 



where 



F{x\T) := -Tin J 



H(x,y) 

dye T . 



This is a standard adiabatic limit. 
Note two more properties. 
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1. If our whole system (protoplasm) is described by a Gibbs distribution: 

1 H(x,y) 

w(x,y) = -e t , (19) 

then, the distribution of probability that protein molecules are situated in the given point of 
phase space can be achieved by integrating ( fl9|) by dy. The distribution of probability w(x) 
to find protein molecules in the given point of the configuration space is: 

F(x\T) 

w(x) = conste T . (20) 

That is, again we received the Gibbs distribution in which the Hamiltonian is the effective 
Hamiltonian for protein molecules we received above. 

2. By using a standard formula we can calculate the free energy of protein system F'(T) 
for Hamiltonian F(x\T). Elementary calculations give: 

F'(T) = F{T) : = -Tin J dxdye'^ , (21) 

i.e. F'(T) equals to the free energy of the whole system F{T). 

Now let's define the form of our effective Hamiltonian F(x\T) as a function of coordinates 
and momentas of protein molecules. We suppose the contribution to F(x\T) nontrivially 
dependent on x is caused by Van der Waals interaction between protein molecules and 
therefore F(x\T) is given by: 

F(x\T) = E km (p) + Fo(T) + J2 V ^ - <Zil T )> (22) 

i>j 

where Eki n {p) is a kinetic energy of protein molecules as material points, Fq(T) is a function 
of temperature and g» pi, i = 1,2,3, ... are Cartesian coordinates of protein molecules and 
canonically conjugated momenta. V(x\T) — is a pair interaction potential of the following 
form: 

V(q\T) = +oo if |g| < 2r , 

V(q\T) = ~, if \q\ > 2r , (23) 

where C is some positive constant. The explicit formula expressing C in terms of e, e', tq can 
be retrieved, for example, from (Lifshitz and Pitaevsky, 1978). Now, we find thermodynamic 
variables for the case just described. The solution of this problem providing that 

min 2ro<r |V(r|T)| 



T 

12 



< 1 (24) 



is stated in many textbooks, for example (Landau, Lifshitz, 1995). However, in real cells (for 
example, erythrocyte) this condition is not fulfilled; this is a subject of a special analysis in 
the next section. For simplicity here we suppose that condition ( l24|) is fulfilled. 

In the case of the living cell, as we shall see later, the free energy of a cell is given by 
expression: 



F(V,T) t = F (V,T) + F id (V,T), (25) 

where Fq(V, T) is a free water energy where all the proteins are eliminated and F^ is free 
protein energy calculated as if it was an ideal gas. For the dead protoplasm 

F(V, T)i = F Q {V, T) + F ld (V, T) + AF(V, T), (26) 

where 

AF(KT) = -THj--J d ^^- i ^ V(q '- , ' IT \ (27) 

V V 

Note that when we use formulas (1251) . ([26"]) . (|27|) . we neglect the conformational part of the 
free energy. There is one more omitted contribution to the free energy, which is caused by 
possibility of protein molecules rotation as a unit. However, these contributions for fully 
unfolded and folded protein conformation differ only by the value of const T and therefore, 
as it will be clear from the following, omitting of those contributions does not influence the 
final results. 

For the case when (EMI) is fulfilled and a gas is so much rarefied that we can take into 
account only pair collisions, the calculations of integral ( 12?]) are made in many textbooks (for 
example, see (Landau and Lifshitz, 1995). But to find AF(V,T) in the real case we should 
derive an expression for AF(V, T) when condition ( 1241) is fulfilled and only pair collisions 
are taken into account. Therefore, we give here derivation for AF(V,T) (condition (I24p is 
fulfilled). Put by definition (U(q\T) := £ V(ft - qj\T)). Then 

l<i<j<N 

AF(V,T) = -Tl»[/.../ ^....f^{e-^ - 1} + 1] (28) 

V V 

If we take into account only pair collisions and suppose them rare, then the whole config- 
uration space of the system C = M. 3 x ... x M. 3 should be divided into subareas with equal 

13 



volume Qj, 1 < i < j < N such that in each of them collisions of i-th and j-th particles 
happens. But in C id U{q\T) = V% - qj\T) So: 

d 3 qi d 3 q N r nn-ym 



mV,T) = -TH £ /^....^{e-^^-D + l]. (29) 

l<i<j<NJ 

But pairs (i, j) such, that 1 < i < j < N could be chosen by N ( N ~ 1 *> « ^ ways. So, we have 



AF(\/,T) = -Tln[— / / d s gi d%{e - 1} + 1] 



v v 



-Thl/^-D + H (30) 
y 

Expanding the logarithm in a Taylor series near 1, the final result is 

AF(V,T) = -T^ J d'qie-^ - 1}. (31) 
v 

In the area r := \q\ < 2tq {e — 1} = — 1 and in the area r > 2tq {e — 1} = 

■ V ^ T ^ due to weakness of interaction. Put by definition: 

3 ' 



a = 2n j \V(r\q)\r 2 dr. (32) 

2r 

Notice by the way that using our choice of potential 

a = 12? (33) 
With these notations we can rewrite F(V, T) in the following way: 

AF(V,T)=T^(b-^). (34) 



So 



N 2 

Sd — 5*0 + Sid — b-^--< 

E d = E + E id - (35) 
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where S d and E d are entropy and energy of the dead protoplasm respectively, S and E are 
entropy and energy of the protoplasm where all the protein molecules are eliminated and S id 
and E i( i are entropy and energy of the ideal gas of protein molecules. 

That is a definition of the thermodynamic features of the dead protoplasm in our model. 
Now, let's consider the living protoplasm. Let's divide the effective Hamiltonian of the 
protein system obtained by the method of adiabatic limit into two summands: 

H = H v _ d _ w + H', (36) 

Here H v _d-w contains the kinetic energy of proteins as material points and energy of the Van 
der Waals interaction between them. The Hamiltonian H' depends on variables describing 
internal degrees of freedom of proteins. We suppose that in a certain sense H' <C H k , 
where is the kinetic energy of all protein molecules. Next, according to our common 
view (Prokhorenko and Matveev, 2011), in the living protoplasm some first integrals in the 
involution Ki, K n are active and a statistical weight of protein molecules is given by: 

/N n 
II dpidqt J] 5(Kj - K^)8(H V ^ V + H' - E), (37) 

i=l j=l 

where Pi, qi are momenta and coordinates of i-th protein. 

Let's try to define the form of integrals K^. According to Ling, intracellular water in the 
resting living protoplasm is in the bound state and protein molecules form a paracrystals. 
Since we suppose the living state differs from non-living by activity of integrals K iy then 
it's reasonable to assume that fixation of protein molecules in points of lattice is performed 

n 

by means of multiplier Yl ~ Kj) i n the integral of the statistical weight definition. 
Therefore, we just suppose that 

n N 

n^-^)=n^-^, (38) 

3=1 i=l 

where are coordinates of i-th. point of protein molecules lattice. But on the support 

N 

of Yl 5(qi — q'j) the potential energy of protein interaction V( x i ~ x j\T) is a constant. 

i=l i>j 

N 

Furthermore, we can suppose that on the support Yl — Qi) the potential energy of 

i=i 

proteins ^2V(qi — qj\T) = 0. Indeed, if proteins are situated in points of the lattice we 

i>j 

mentioned above, then 

J>(%-^)~^, (39) 

i>j 
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where r is a minimal distance between proteins. But r ~ (jf) 1 ^ 3 - Therefore: 



J2v( gi - qj \T)~N(y) 2 . (40) 

i>j 

Further we show that if the potential energy of protein molecules interaction is neglected, 
then, when the cell is dying, the quantity of emitted heat and the number of potassium ions 
released from the cell, calculated for one protein molecule, is a linear polynomial in y with 
an accuracy up to logarithmic factors. So, as follows from ( )40|) . in ( l37l) in the limit of low 
density the potential energy of protein interaction in H v _d-w can be neglected just kinetic 
energy H^. Eventually, for the statistical weight: 

r N N 

W(E)= l[dp l dq i l[5(q J -q' :j )5(H k + H' - E). (41) 

J i=l j=l 

Since H' <C we can write over this expression in the following way: 

/N N N 2 

J] dPidqi J] - - E). (42) 

i=l j=l i=l 

Consequently with an accuracy of inessential multiplier, the statistical weight W equals to 
the statistical weight for the ideal gas. Therefore, in our model thermodynamic variables for 
the living protoplasm take the form of: 

Ei = E + E id , 
Si = S + S td . (43) 

Here E\ and Si are entropy and energy of the living protoplasm. Obviously E& < Ei, i.e. 
our model predicts that activation and death of the protoplasm is an exothermal reactions. 
Numerical evaluations are given in the next section. So, we have the following amount for 
the quantity of released heat: 

Q - ^- (44) 

The ability of variables qi or closed to them play the role of motion integrals in our model 
is appeared from the following additional conclusions. Let's denote by Ek in a kinetic energy 
of one protein molecule. Obviously 
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Therefore, in the limit M — > oo (very large mass of protein) q { = and are motion 
integrals. 

Now we define the number of potassium ions releasing from the dying cell. In our previous 
work (Prokhorenko and Matveev, 2011) the following formula was obtained. Let's suggest 
that the number of first integrals in the involution is so large that the number of active 
integrals can be characterized by continuous parameter s G [0, 1]. In addition, the number 
of active first integrals is increasing function of s and a case of s = corresponds to the 
case when none of integrals is active, and case of s = 1 corresponds to the case when all 
the first integrals are active. Supposes s infinitesimally varies s i— > s' = s — 5s, 5s > 0, 
where s is infinitely small. Now, let's define the function of entropy (5f)(S) by the condition 
(Sf)(S) = (SS)e, ( s h-» s' = s — 5s). This definition is correct due to a proven assertion 
(Prokhorenko and Matveev, 2011) that (5S) E ( s h-> s' = s — 5s, 5s > 0) is a constant along 
any adiabatic process. So, let s i— > s' = s — 5s, 5s > 0, 5s, is infinitely small. As it was 
shown in (Prokhorenko and Matveev, 2011), the increasing of number of potassium ions in 
the cell 5N can be found using the following formula: 

5N = T(5f)'(S)(—) T , (46) 

where /i is a chemical potential of potassium ions in the cell. 

First, let's define (§^)t- Suppose by definition F (V,T) is the free energy of the proto- 
plasm where potassium ions absent. 

Suppose by definition F (V, T) is the free energy of the ideal gas at temperature T while 
its mass equals to the mass of potassium ion M and the gas contains only one molecule. It 
is known (Landau and Lifshitz, 1995) that 

MT 3 

<t ) (V,T) = -T\n{V(—y>]. (47) 

Suppose by definition ip(V,T) is potassium solvability in the protoplasm, i.e. change of the 
free protoplasm energy during transition of one potassium ion from the infinity to the given 
point. Then the change of the free protoplasm energy after adding one potassium ion to this 
protoplasm is just a sum of (p(T, V) and ip(V, T). So, the free energy of all the protoplasm is 

F (V, T) h. F(V, T) = F (V, T) + <p(T, V) + ^(V, T). (48) 

If there are N potassium ions inserted to the protoplasm, but iV is still very small, then an 
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interaction between different potassium ions can be neglected and we have: 

- N 
F (V,T) ^F(V,T) =F (V,T) + N<p(T,V) + NiP(V,T) + NT\n(—). (49) 

e 

The last summand arises because of a common combinatorial multiplier included to 
definition of the partition function of potassium ions. 

We neglect the solvability of potassium ions in our calculations, in other words we accept 
1>(V,T) = 0. 

Put by definition 

Eventually, after introduction of N potassium ions into the protoplasm, the free proto- 
plasm energy changes as follows: 

F (V, T) H- F(V, T) = F (V, T) + A^(T, V) + NT ln( ^ ). (51) 

Notice that the quantity A( ^ T ^ <C 1 since by order of magnitude this value is a density of 
particles in phase-space cells of constant-energy surface and thus it is very small, because 
the gas of potassium ions, as we admitted, is strongly rarefied (particularly, this allows us 
to use formulas of the classic statistical mechanics). 

Potassium adsorption is taken into account in the solvability ip(V,T). We supposed that 
this solvability equals to zero but this assumption is not quite right. After cell death the 
chemical potential of potassium ions, obtained by differentiating the free energy over with 
number of potassium ions, remains constant. Since while the cell is dying, it issues ~ 0.98 
of the total amount of potassium ions, solvability ip(V,T) increases by Tin 50. Though this 
section contains exact equalities, actually this equalities are approximate because we neglect 
a unity compared to ln( e ^^ r J and change of ip(V,T) by value of Tin 50. However, the 
following analysis establishes that contribution of those omitted summands is inessential. 

When A" potassium ions are inserted into the protoplasm,the entropy changes as follows: 

S Q (V,T) H- S(V,T) = S Q (V,T) - AHn ^^ ). (52) 



We have 



( dS (||)r , v 
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Further: 



And: 



Therefore: 



Eventually: 



and 



W^-^W^- (54) 



,dF(V. T, N) , , , iV x 



,dS. N. ,eA(V,T), 



< yjv = jv(W(g)in[ &A( ^ r) ]. (5? 



4 Van der Waals Model of Protoplasm. Numerical 
Evaluations 1. Potassium Ions Efflux from the Cell 
and Heat Release 

This section is concerned with getting numerical evaluations basing on the Van der Waals 
model (see above) and comparing them with experimental data. To specify model parameters 
(cell size, quantity of proteins and ions in the cell) we chose a human erythrocyte, the well- 
studied cell having a relatively simple structure-function organization. 

It's known that potassium ions density in the living erythrocyte is estimated as n = 
6.02 x 10 19 cm -3 . Let's assume that when an erythrocyte is dying the potassium concentration 
in this erythrocyte becomes equal to the potassium concentration in the blood plasma, 
2 — 4mmol/l, i.e. about 0.96 — 0.98 of potassium ions, which the living cell contained, 
release from the erythrocyte. 

We consider the cell at a temperature of 300JT. A potassium nucleus contains 19 protons 
and 20 neutrons, therefore, the mass of potassium ions can be estimated as 39m„ where m n 
is the mass of neutron. 
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In the previous section the formula was derived describing the change of a number of 
potassium ions in the cell when the infinitesimal change of a number of active integrals 
s i— > s' = s — 5s, Ss is infinitely small. Let's write it one more time: 

8N = N(8f)'(S)H eA ^ T) }. (59) 

The last formula ( 158]) can be interpreted as a differential equation on N. Being integrated 
this equation expresses of potassium efflux from the dying cell through other cell parameters. 
First, let's find (Sf)'(S). We have: 

(99 

Sf(S) = (SS) E = (SS) T - (5E) T (—) T . (60) 



But 



,dS_ 1 
> 9E )n t 



(61) 



Therefore 



T 

Using the fact that (see Prokhorenko and Matveev, 2011 

d{5S) T 1 



5f{S) = (5S) E -h5E) T . (62) 



d(5E) T T' 



(63) 



we find 



d5f(S(T)) 1 

dT = Y^ 5E)T - (64) 

And finally 

WW = y^ 6E)t % = ^t {6E)t ' (65) 

where Cy is the heat capacity of the cell at the constant volume. Now, we assume the heat 
capacities of living and dead cells are almost the same, and in all following formulas the heat 
capacity of a cell can be replaced by an average value Cy which does not depend on the 
number of active first integrals in the involution. Eventually, we have 

f = -m {SEilH ^) v (66) 
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For convenience of further calculations let's introduce a new variable 



A:= 7A^ry (67) 



Then we have: 



51nA=-^i-(lnA)(5£) T . (68) 
As before, we use lower indexes I and d to denote variables relating to the living and 
dead cell respectively. Equation ( 1681) is easy to integrate resulting: 

= -P^), (69) 
In A ; GyT 

where Q is a heat amount released from the cell while it is dying. 

The last formula gives an implicit expression for heat generation Q and its derivation was 
based just on common thermodynamic considerations without regard to properties of any 
certain model. Therefore, it can be used for verification of our generalized thermodynamics. 

However, note that (169 j) is inconvenient for calculations, so, let's simplify it using small- 
ness of A;. For this purpose let's consider the expression In x, where a; is a very small positive 
number. Then In a; value is very large in modulus and has a sign minus. Let's increase x 
by factor of k where k is not too large natural number lnx i— y In a; + Ink i.e. practically 
unalters. We have: 

I lnA d | 

exp{ / -}. (70) 
J x 



In A rf ,lnA„ 



In A ; In A/ 

I In Ail 

According to the newly stated remark, on the whole integrating interval we can replace x 

by | In A; | and eventually we obtain: 

In Ad , A; , i , , 

InA, l A/ 1 J 

Eventually, we receive the following implicit expression for heat generation: 

^ex P {|.nA,|JL } . (72) 

Data listed in this section are enough to calculate | In Aj | . Omitting corresponding numerical 
calculations, we present the result: | In A; | ~ 16.1. The resulting formula for heat generation 



is: 



W d = eXp{ C^T } - (73) 
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From here it is easy to find 



C V T- N d 



j± fa 50. In 50 = 3.93. By taking a logarithm of both parts 



we find: 



Q 



0.24. 



(74) 



C V T 



To check this heat generation value we can carry the following qualitative reasoning. Let's 
consider a human cell at temperature Th = 37C = 310K or 37°C. The room temperature 
(normal for a man's functioning) T r = 20C = 293K. According to Ling's theory, the cell 
life activity is expressed as the cycle motion " resting -B-exciting" . If the cell was always in 
thermodynamic equilibrium with the thermostat, then, at transition "resting^-exciting" , the 
Q heat would released. Let's accept that in the resting state a cell has a room temperature. 
Let a be a digit making Q = aCyT '. If the cell is heat sealed, then in transition from the 
resting state to the excited state the temperature of the cell rises from T r to T e and it is 
likely to be correct that there is an approximate correlation (T e — T r )/T h = a. Now, as the 
cell moves through time cyclically, we have a mean temperature for time T h . It is reasonable 
to suggest that = (T e + T r )/2, or in other words T e — T r = 2(7^ — T r ). But T^ — T r = 17K. 
Therefore T e — T r = 3AK and a = (T e — T r )/Th = 34/310 = 0.11, which coincides with our 
result in the order of magnitude. The fact that the calculation for heat emission of the dying 
erythrocyte using data on potassium efflux from the cell appears approximately 2.5 times 
higher than just stated value is reasonable because the exciting can be considered as a stage 
towards death (Matveev, 2005). 

The assumption we stated that the erythrocyte temperature at rest equals to the room 
be explained in the following way. Since the erythrocyte moves through time 
cyclically: restingoexciting, then the erythrocyte can be considered as a heat engine, and 
room temperature is a temperature of a cooler of this engine. As the room temperature is the 
most comfortable for a man, we may consider that the erythrocyte at the room temperature, 
as at the cooler temperature, operates in the most optimum mode. The resting temperature 
T r is the lowest temperature achieved by the erythrocyte during all the resting -B- exciting 
cycle. If T r > T room is correct, then heat transfer from the erythrocyte to the cooler will 
beat tended by heat transfer from a warmer body to a colder one, i.e. entropy increase. This 
means the erythrocyte as a heat engine would operate not optimally. Conversely suppose 
that T r < T room . If the erythrocyte is functioning in the optimum mode (without entropy 
increase), then the erythrocyte have to pass the part of the cycle when the erythrocyte 
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temperature T <T r being surrounded by an adiabatic "cover" . It becomes incomprehensible 
why does the erythrocyte need this part of the cycle and also, according to Carnot's theorem 
on the efficiency of heat engines, erythrocyte efficiency could be raised by means of the cooler 
temperature (room temperature) decreasing. 

Now let's try to calculate Q basing on our investigated (Van der Waals) model. For this 
purpose we should find an interaction constant C in the law V(r) = -fi. To do this, we 
proceed from the formula taken from (Lifshitz and Pitaevsky, 1978). Suppose there are two 
parallel non-overlapping semi-spaces and a distance between them is /. Let e be a dielectric 
permittivity of the semi-spaces and e' is a dielectric permittivity of the cavity between them. 
Suppose I is such a large that if w is a frequency of the electromagnetic wave distributed in 
the cavity between semi-spaces and having wavelength /, then < T. That is obviously 
our case. Then, the force of attraction between two semi-spaces as per unit area of the 
border space of each semi-space is equal to: 

+00 

p =^S^& e '-^ ix - (75 » 



Now, suppose that the semi-spaces consist of folded hemoglobin molecules and the cavity 
between them is filled with water. The dielectric permittivity of the water e — 81 and the 
dielectric permittivity of the hemoglobin e' ~ 2. Therefore, the factor {f^?} 2 ~ 1 and the 
force of attraction between two semi-spaces as per unit area of the border space of each of 
them is equal to 

+00 

P= lLp j A^-ir'dx. (76) 


The involved integral can be easily calculated with any desired degree of precision, for ex- 
ample: 

+00 +00 

-!]-'*= /^[i -«-*!-'* = 



+00 

= J x 2 [e~ x + e' 2x + e~ 3x + ]dx = 



= 2[1 + ^ + ^ + ....] = 2[1 + 1/8 + 1/27 + 1/64 + ...] « 2.35. (77) 
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Therefore, we come to the following formula for the force of attraction between the semi- 
spaces: 

P - 2 ' 35T (7R) 

P ~ To^p- (78) 

Let Vih be a volume of one hemoglobin molecule. Now, we can calculate a potential energy 
U(l) per unit of the surface plane using a formula for attraction potential of two hemoglobin 
molecules ff2"3"j) in the following way: 

+00 +00 +00 

Um—^Jr*]*]*, 1 (79) 



Omitting rather trivial integrating we find 

= -ise 4 (80 > 

Differentiating the last expression with respect to I we find: 

vrC 1 

x 



6V? H P 



Comparing the last formula with (178!) gives the following result: 



C = V? H 7 -^r. (82) 

It follows 

7 05 

a = V 1H T x —— = 9.79 x 10~ 2 V 1H T (83) 

Note that everywhere above we measured the temperature T in energy units. Let T' 
be an absolute temperature expressed in Kelvin degrees. There is a relation T = ksT' 
where ks = 1.38 x lO^^ergK' 1 . Hereafter Q v means the heat released from a cell within 
the framework of the Van der Waals model. We want to estimate Then, for heat 

generation: 

O v = 1.35 x lO- 17 N 2 &{^)erg x K~ X T' (84) 

V V H 

where Vh is the volume of all hemoglobin contained in the dead erythrocyte. 

We need to know the cell volume V, volume of all hemoglobin Vh contained in the 
dead cell, mass of all the hemoglobin Mh, specific heat capacity of hemoglobin per unit 
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mass C M - Calculating this data we find V = 10~ 10 cm 3 (Levine et al., 2001). C M = 3.2 x 
10 7 erg g^K^ 1 (Kholodny et al., 1987). Further it's known that hemoglobin concentration 
C H = 5mmol x L" 1 , C H = 5 x 10~ 6 moZ x cm" 3 (Van Beekvelt et al, 2001). The volume of 
"dead" hemoglobin is V H = 32.6 x 10 3 cm 3 /mo/ (Arosio et al., 2002). N = 3.02 x 10 8 (Van 
Beekvelt et al., 2001). Therefore, 

= 0.163 (85) 

But hemoglobin concentration in the human erythrocyte is 0.33g/cm 3 (Van Beekvelt et 
al., 2001) and its density in the dead cell (in the supercluster) is approximately two times 
greater than the water density (Van Beekvelt et al., 2001; Arosio et al., 2002). It follows 
the mass of the erythrocyte M c = 1.16 x 10~ 10 g. The heat capacity Cy is Cy = M c Cm = 
3.71 x 10~ 3 erg K~ x . 

Let's calculate N 2 ^f. Taking into account that N 2 = 9.12 x 10 16 , we find 



N 2 V, 



H = 1.49 x 10 16 . (86) 



V 

Eventually 



Qv 54^ (87) 



CyT V H 

Here we see that if we take a true value 3.3 x 10~ 9 for then the resulting value for 
-Q^f, is many times smaller than = 0.24 we found before. 

This fact has the following explanation. Our derivation of corrections to the free energy 
was correct only upon the following condition: 

min 2ru<r ^(r|T)| ^ y (gg) 

Let's calculate the value in the right part of this inequality. Suppose two protein molecules 
are situated so that a distance between their centers is slightly more than 2r . Then, the 
potential energy between them is U — — C ^ 2? ^ 6 . We have shown above that C = V 2 H T^§-. 
Here, after elementary calculations, we find 

^ = ™5 w 2 . 5 x 10 -2 (89) 
T 288 ' v ; 

i.e. our criterion is really fulfilled. However, when two protein molecules come to each other 
such close that 2nhj 1, the forces whose contribution we did not take into consideration 
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before begin to play a significant role. We mean Casimir forces caused by the energy of 
zero-point oscillations of the electromagnetic field in the space between proteins, and which 
significantly exceed forces we took into consideration before. These forces can be interpreted 
as chemical ones. As a result, considering these new forces, it is reasonable to expect that in 
the real dead protoplasm the protein molecules stick together to balls or superclusters, and 
we should consider this effect in derivation of AF(V,T). The following section generalizes 
the derivation of corrections to the free energy in case of possible clustering, and it becomes 
clear that if V KH denotes the volume of such a cluster and Q v is a heat generation of the Van 
der Waals model, and Qkv is a heat generation for the same model (considering the possible 
clustering), then 



As a result, we have 



Qkv = Qv^-- (90) 



Q kv m Vkh (gi) 



C V T V H 

If this new formula considering collective phenomena is correct, then we shall find that 
the dead protoplasm should contain ps 225 superclusters or aggregates, the properties of 
which we cannot characterize yet because the conclusion about the number of clusters is a 
result of a rather common theoretical analysis. On the other hand, the protein aggregation 
in case of the cell death is a well-known phenomenon. 

5 Van der Waals Model of Protoplasm. Numerical 
Evaluations 2. Clustering of Protein Molecules 

In this section we estimate a correction to the free protein energy AF with respect to their 
clustering. The correction AF(V, T, N) is 

AF = -Tln|^i...^e-^. (92) 

Here ?i,...,?jv are Cartesian coordinates of all molecules having numbers 1,...,N and 
symbol q means coordinates of all molecules: q := (gi, ...,q N ). U(q) is given by 

U(q) = U u (q) + U s (q). (93) 
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Here U u is a potential energy of the Van der Waals attraction between particles: 

l<i<j<N 

(94) 

where 

V u (q) = -Cp. (95) 

Hereafter C value is considered as a small parameter used to make any asymptotical expan- 
sions. 

U s {q) = U 1 s {q) + U 2 M ) (96) 

where U] (q) is a potential energy of repulsion between protein molecules, arising due to their 
volume is low bounded. For example: 

l<i<j<N 

V s l (q) = 0, if|g|>2r , 
K 1 (g) = +ooif|g|<2r , (97) 

r is a radius of the protein molecule. U% is a potential energy of Casimir forces, arising in 
very closed distances between protein molecules due to zero-point oscillations of the electro- 
magnetic field in the gap between separate proteins. These forces are short-ranged but they 
have a high degree of cooperativity which contributes to clusters formation. 

Potentials U s , U], U 2 S we sometimes call superpotential, and forces corresponding to 
them, superforces, as U}, U 2 are much greater than U u (q) and the effective consideration 
of U s in the Gibbs exponent comes down to the fact that the whole configuration space of 
molecules is replaced by its part. 

Let's represent AF in the following form 

AF = AF + AFi + (98) 

where AFq is a zero-order value of vanishing on C, AF± is a value of the first order of 
vanishing on C and so on. 

Let's begin from the definition of AF , in other words, put C = and U(q) = U s (q). 
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Due to presence in U s all the molecules stick together to clusters (balls B) of N B molecules 
in each. Let N k be a number of such balls (clusters). Obviously N = N B N k . 
By definition all the clusters identical and there are 

I N\ 

(99) 

of equivalent ways to distribute molecules by clusters. Let's fix one of these ways which we 
will follow hereafter, where molecules having coordinates qi, qN B belong to the first cluster 
Bi, molecules having coordinates qN B +i, q2N B belong to the second cluster B 2 and so on. 
During F calculation we naturally neglected interaction between clusters, so: 

N\ f 1 f d 3 q x d 3 q A r E ^ u b (ii,^,in b ) 



AF = -Tin— -[— - / — ^H^, (100) 

N k \ [ N B \ J V V 1 ' K J 

where by definition we set: 

u B ( gi ,...,q NB ) = v s(Qi-Qj) + Us(qi,-,qN B ). (ioi) 

l<i<j<N B 

We proceed from the assumption that the cluster is a ball and bonding forces between 
molecules in the cluster are so strong that the cluster volume equals to a sum of volumes of 
molecules it consists of (that why cluster density can exceed the water density). 

Further, one of variables q±,....,qN B is a center of inertia of cluster B. The choice of such 
a variable can be performed in N B ways. Supposing the center of the cluster inertia is q± we 
find by integrating it: 

TV 1 f d 3 a-> d 3 aM u B (o,..., qN ) 
AF = -Tln^-[ 7 —— / ^....^* e ^^} Nk , (102) 

And the center of cluster B is 0. 

Then,we assume every molecule in the cluster moves in an effective field W and 
Ug (0, . . . , q 7v s ) = W(N B — 1). Instead of integral J d 3 q 2 ---d 3 qN B the presence of mul- 

Bx...xB 

tiplier e t imposes the use of integral 



d 3 q 2 ...d 3 q NB , (103) 

Bx...xB 

where prime points at the incompressibility of molecules. Using AH (^r) N we receive: 

1 r 1 ,N Nb f d 3 ^ d 3 q NB ^-w(N B -i)^ Nk 



AF °- Tla m l (N^Ty<- ) a J ^--v^t-^r. (104) 

Bx...xB 
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But integral J d 3 q 2 ..-.d 3 q NB is easy and equals to 

Bx...xB 

i 

J d%....d 3 q NB = (N B - m^) NB ~\ (105) 
Bx...xB 

where Vb is the cluster volume and := is the volume of one protein molecule. As a 
result, we obtain: 

Vr \t i 1 -W(Nr-I) 

AF = -TiV fc lniV i? [(^)^- 1 ^- T e^^], (106) 

where V& := j^- is the volume which falls on one cluster. 

Now, let's derive the equilibrium condition for the cluster B from which we essentially 
can define the volume of this cluster. Let's consider one equilibrium cluster in the volume 
Vk- Now we are generally interesting in configuration part of the free energy since it splits 
from the kinetic one. 

Let's add one more protein molecule to the volume Vk- Its free energy is 

A = -Tln j y/q = ®- (107) 

If this protein falls to the cluster, then Nb is increased by a unit, but since the cluster is 
equilibrium AF should not change. In other words, the following equality should take place: 

^ = °- (108) 
dN B 

The equilibrium condition fl 1 8 [) is actually a common thermodynamic condition of equi- 
librium if we suppose a number of protein molecules in a cell N to be a parameter which 
has an influence on the system condition. But we can use this condition only if iV is not 
preserved, put it otherwise, is not a motion integral. But the fact that iV is not a motion 
integrals follows from the fact that in the living cell the quantity of protein molecules does 
not remain invariant because of permanently proceeding processes of proteins synthesis and 
degradation. 

The equation (1 108[) leads to the equality: 

I V R V h W 

— + ln(-^) + (N B - 1)-^ - 1 - - = 0. (109) 
i\b Vk Vb i 
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In other words: 

- = er. (110) 

Finally: 

AF = -TN k \nN B + T[N - N k }. (Ill) 

As the number of clusters is small (N k <C JV), then in square bracket of the second summand 
N k can be neglected as compared with N. Since Nb conversely is very large, we can neglect 
a logarithmic term as compared with TN . As a result, we find: 

AF = TN. (112) 

This term just comes down to the additive renormalization of the entropy and therefore in 
can be omitted. 

Now, let's calculate AF 1 (V, T, N). We have: 

AFi = AF — AF + 0(C 2 ). (113) 

Let's recall, that C is considered as as mall parameter. It won't be a great mistake if 
instead of AF we take AF calculated by iV — 2 molecules in the same volume V at the 
same temperature T. We have: 

AFi = -Tin J <^.J^L e =m±^ + oiC 2 } 
= -Tln/^....^ J] [e-^^-l + l^^^+O^ 2 ) 
= -Tln/^....^{l+ J2 [e~ V -^-l]}e^^ + 0{C*) 



l<z<j<N 

= -T H l + [ ^....^{ E - 1|e ^r^ } + o(&). (114) 



l<i<j<N 



N 2 f d 3 qi d 3 qN r Vu(n-q2) i1 _ -^w+ AF n - 1 , ,, , , » 



Eventually, we obtain: 

AF 1 = -Tln{l + — y r_^....r_^ [e — ^ _ i ]e -^— } + 0(C a ). (115) 

Since the number of clusters N k is still sufficiently large, we neglect cases when qi and q 2 are 
in the same cluster and suppose them lying in different clusters Bi and B 2 . Values qi and 
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q 2 are simultaneous centers of the corresponding clusters with probability -b-. Taking this 
into account: 

Afl = _ rln{1 + ^|^....^ K e-^)-l]e^ i } + (^). (116) 

Here variables q\ and q 2 are already centers of inertia of clusters B\ and B 2 . Variables q[, q 2 
are equiprobably and independently distributed along clusters Bi and B 2 and angle brackets 
mean the averaging over positions q[ and q' 2 . 

Considering the meaning of AFq, the last equation can be rewritten in the following way: 

AF, = -Tl„{l + ^/Vl / rfVrffc-^) - 1] 



N Nb ~ 1 f d 3 q 2 d 3 q Ng -t/ a (o,-, TO ) l2 2 



|<?|>2r s 



L (N B -1)\ J V V 

Bx...xB 

AF l = -T\n{l + —N 2 B J d 3 q[(e ^} - 1] 

kl>2r s 

x[(^- 1 e-V ! - i ]2} + (C7»), (117) 

where r B is a radius of a cluster V B = ^fr B . But taking into account the clusters equilibrium 
conditions we derived, the expression inside the last square bracket equals to 1. Eventually: 

AF 1 = -Tln{l + ^iVl J rf 3 g[(e- M ^)-l]} + 0(C 2 ) 

kl>2rs 

= - T W N ° J ^[(e-^>-l]+0(C 2 ). (118) 

|?|>2r s 

This formula assumes that cluster B\ has a center of inertia in zero, q is a coordinate of 
center of inertia of cluster B 2 , q[ and q' 2 are variables which are uniformly and independently 
distributed over cluster B\ and B 2 respectively. 

This formula particularly proves the rule (|90|) quoted in the end of the previous section. 

Now, let's define a part which "dead" protein in the protoplasm occupies in the whole 
cell volume. Suppose the cluster of dead protoplasm proteins is situated in the volume and 
has an V := jj- equilibrium volume. If the cluster has the equilibrium size, then inserting 
one protein molecule to it needs a zero work A, in other words, an average potential energy 
of the protein molecule in the cluster equals to zero. But all the protein molecules are 
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situated in the neutral external field outside the cluster and in the field W inside the cluster. 
Let's consider the chosen protein 25, it has three degrees of freedom described by Cartesian 
coordinates q 1 , q 2 , q 3 and radius q A (the protein can be compressed). The only thing we can 
do is to take a positively defined quadratic form of coordinates q 1 , q A as a raw estimation 
for the potential energy of protein inside the cluster V^q 1 , q A ): 

V(q 1 ,...,q A ) = W+ 4), (119) 

i,j=l,...,4 

where qQ,..-,q 4 are equilibrium coordinates of the protein molecule. For kinetic energy of 
the chosen protein 03 as a raw estimation too we take a positively defined form of velocities 

•1 -4 

q ,...,q : 

E km = J2 ( 12 °) 

ij=l,...,4 

Note, that protein molecule in the water does not behave as an oscillator along g 4 because 
a water is incompressible. 

But in the case of classical mechanics the average potential energy of one-dimensional 
oscillator equals to |T. Therefore, the condition A = leads to W = —IT . This equilibrium 
condition cluster lets us find: 

^ = e f = e -2^ 0.136, (121) 
which is in a good agreement with true value ^ = 0.163 (see section 4). 

6 Superfluid Bose Gas on Protein Configuration Space 
as Model of Living Cell Protoplasm 

In this section we discuss the second model of the living cell protoplasm, the superfluid Bose 
gas on the protein configuration space. The common description of this model is given in 
the introduction section. While investigating this model we'll show that when the Ling's 
protoplasm dies, the folding of protein molecules occurs, how the Ling's theory postulates 
on the qualitative level (Ling, 2001). 

The protoplasm is considered at zero temperature. This allows us to apply essential 
tools of theoretical physics but conclusions made for such extreme conditions are useful for 
understanding the properties of the real cell. 
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Let's proceed to conclude our "superfluid" model. Suppose there are N protein molecules 
in the protoplasm, x iy i = 1,...,N are coordinates of their centers of inertia and <jj are 
parameters, passing some space X with measured d\x describing inner freedom degrees of the 
protein molecule (its configuration). Then, the state of N protein molecules is described by 
the wave function ^(x±, a±, x^, cr N ) depending on coordinates x iy Oi of protein molecules. 
This wave function fulfills the normalizing condition: 

J \V(x 1 ,a 1 ,...,x N ,a N )\ 2 dx 1 ,...,dx n dfx(cr 1 )...d(x(<T N ) = l (122) 

and symmetry condition VP G Sjy permutation of the set of N elements 

(P^)(x 1 , <7i, ...,x N , cr N ) : = t(x P (i), <7 P (i), x P ( N ), a P ( N )) = x(P)^(xi, -, %n, <t n ),(123) 

where x(P) = 1 if protein molecules satisfy to the Bose - - Einstein statistics, and 
x(P) — sg n (P) (signum of permutation) if protein molecules satisfy to the the Fermi - 
Dirac statistics. 

To describe this system, by using common considerations we'll try to find its Hamiltonian 
as per the adiabatic limit method (see the previous section). Suppose V(xi, <Ji, xn, ctjv) 
is a change of the smallest eigenvalue of the protoplasm Hamiltonian while adding there iV 
protein molecules having coordinates x±,ai, ...,xn,ctn, and V is an operator of multiplica- 
tion of wave function of our system on V(xi, a±, xn, <Jn)- Then, obviously, the effective 
Hamiltonian of our system is 

N 

H eff = J2T t + V, (124) 
i=i 

where is a kinetic energy operator of i-ih protein as a material point: 

fi = —Vl (125) 

2M 1 y ! 

To simplify our model let's suppose V describes only a pair interaction of protein molecules. 

N 

^ = E^+ E Vij, (126) 

i=l N>i>j>0 

where in Dirac notation 

(a 1 ,x 1 ,...,x n ,a n \V i \x 1 ,a' 1 ,...,x n a' n ) = U^Xi, a^, a-) (127) 
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for some function U^x'^a'^x'^a^) fulfilling obvious conditions arising from the Hermitian 
nature of V$. Similarly 

((Ti,xi, ...,x n ,(T n \V iij \x 1 ,(T' 1 , ...,x n a' n ) = U 2 (x i ,a i ,x j ,a j \x' i ,a' i ,x' j ,a / j ) (128) 

for some function U2(xi,ai,Xj,aj\x' i ,a' i ,x'j,a'j) fulfilling evident conditions arising from the 
Hermitian nature of Vij. 

We consider the protoplasm at the zero temperature (measured in the energy scale). 
The following objection may arise against this acceptance: at zero temperature all the 
metabolic process in cell are stopped and the cell dies (it is the way how the non-equilibrium 
thermodynamics understands the life and the death). However, it is well-known that cells 
(even embryos including human ones) come back to life after freezing in the liquid nitrogen 
when there are no flows of matter or energy at all. Our approach differs: we consider the 
living state as stationary and thermodynamically sustainable, but not the equilibrium one. 
Precisely this state is "freezing" and can keep its properties even at the absolute zero. That's 
why the frozen cell comes back to life when temperature conditions became normal again. 
The thermodynamic identity of the resting state of the living and cells allows us to apply 
the same analytical apparatus. 

As the average kinetic energy of a heat motion of the protein molecule Tj at temperature 
T is equal to kT, so at the zero temperature kinetic energy equals to zero. Therefore, we can 
omit the term corresponding to kinetic energy from the Hamiltonian H e ff. In other words, 
when T = we have the following expression for i? e // : 

N 

= + E Vij. (129) 

1=1 N>i>j>0 

With regard to the protein molecules themselves, according to Ling's model of the cell 
(Ling, 2001), they are situated in points of a lattice and emerged into the volume of water 
associated with these proteins (the water is as well structured and ordered), while i-th point 
has coordinates Xi. ^(xi, a±, ...,xn, un) has the following form: 

N 

o-i, ...,x N , a N ) = C ^ Yl S ( x 'i ~ x P(i) )A p (a 1 , a N )x(P), (130) 
p e s N i=i 

where A p fulfills the following property arising from the property of ^(xi, a±, xn, ctn) to 
fulfill the Bose-Einstein or Fermi-Dirac statistics: 

A pq_1 ((7q(i), ....,<7q(jv)) = A p (a 1 , ...,a N ), (131) 
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and a divergent constant C is chosen from the condition 

r N 

C dx u dx N ([[ S(xi - x-)) 2 = 1. (132) 
i=i 

At T = the whole protein system is in the ground state ^oO^i, oi, x^, &n) which 
can be defined based upon the requirement of minimality of the averaged Hamiltonian H e ff 
over all normalized states. In other words, in the state \l/o the averaged value (^f\H e ff\^f) of 
the effective Hamiltonian reaches its minimum value 

E = mm^\H eff \^). (133) 

Taking into account the aforesaid, in searching the minimum value of (^f\H e ff\^f) we can 
restrict ourselves to states of the form f)130p . Since for states of the form (11301) the coordinates 
of inertia centers of proteins are localized, we may say that acts only in L 2 (X, dfi) and tr- 
acts in L 2 (X 2 , dpi x dfi). Then, our variational problem comes down to finding the functional 
minimum 

f N 
dfx(a 1 )...d l x(a n )A*(a 1 ,...,a N ){(J2v i + }2 V^A}^, a N ) (134) 

J i=l N>i>j>0 

upon the additional condition 



|A((Ti, ...,a n )| dn((T 1 )...dfi((T N ) = l. (135) 
Let's introduce a function V 2 ((Ti, cr 2 [ cr^ , er^a^, x'a) with a formula 

V 2 (ai,a 2 \(r' 1 ,a' 2 \x' i ,x' j ) = U 2 (x' i1 a u x' p a 2 \x[, a[, x' j} a' 2 ). (136) 
Then, in obvious notations 

(a u ...,a N \Vi tj \a' x ,....,a' N ) =V 2 {a i ,a j \a' i ,a' j \x u x j ). (137) 
Let P be any permutation from Sn- Let's define the operator with the following formula 



(<7i, ...,a N \V itj \a' u ...,cr' N ) = V 2 {a h Oj\o\, Oj\x P @, xpyj). (138) 

Let by definition 



= At y V p . 

13 \S N ^ ^ 



VS := > ' V?;. (139) 



p&s n 
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Let's renumber protein molecules in such a way that the protein which starts the numer- 
ation is situated sufficiently far from the cell surface. Then 



1 N 

(a 1 ,...,a n \V i e j \a' 1 ,...,a' n ) = — ^V 2 {a i ,a j \a' i ,a' j \x Q ,x k )] (140) 



N 

k=l 



Let by definition 



1 N 

£/ 2 (<7i,<7 2 K,<72) := — ^2 V 2 (a 1 ,a 2 \a' 1 ,a' 2 \x ,x k ). (141) 



N 

k=l 



Now, we'll show that within the framework of some approximation the class of functions, 
where solution of variational problem we just set up is searched, can be narrowed down to 
symmetrical functions, i.e. such functions that VP G Sn 

A p (a u a N ) = A id {a u a N ) =: A(a ± , a N ). (142) 

We suppose that V it j is proportional to real constant A which is a small parameter. Therefore, 
if parameter A is enough small, the different protein molecules can be considered as non- 
correlated ones, in other words, behaving independently. This means that the wave function 
A((Ti, <Tat) can be represented in the following form: 

N 

A(a 1 ,...,a N ) = H^i) (143) 

i=i 

for a function £(cr) G L 2 (X,dfi) fulfilling the following normalizing condition: 
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Z(a)\ 2 dii(a) = 1. (144) 



Let U l be an operator in L 2 (X,d/i) which have matrix elements in a-representation 
defined by function Ui, and U 2 be an operator in L 2 [X x X, d/i x d/i) which have matrix 
elements in a-representation defined by function U 2 . Let's denote U 2 (£) as an operator in 
L 2 (X,d[i) defined from the following relation: 

(/l^(0b> = </®ei^b®0- (145) 

The approximation we've just described is called a mean field approximation and time evo- 
lution of function £ in it is defined by the following equation: 

i^^V^ + NV 2 ^. (146) 
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Any how, in our approximation 

N 

Aia^.^a^-H^). (147) 

i=i 

As the class of test functions is narrowed down to the symmetrical ones in our variational 
problem, then, instead of ti e ff minimum eigenvalue calculation, we can seek the minimum 
eigenvalue of the Hamiltonian: 

N 

£'e// = £v*+ E Hr ( 148 ) 

i=l N>i>j>0 

Within the framework of our approximation the energy of the ground state of H' e ff is 
equal to the energy of the ground state of H e ff- 

Proposition. There exist functions U'^a^a'-^), U 2 (ai, a 2 \a[, a' 2 ) such as if V'i, V-j oper- 
ators defined by the following relations 

(<7i, ....,a N \V-\a[, ...,cr' N ) = J7(((Tj|cr-), 
(a 1 ,....,cr N \Vr\a' 1 ,...,a' N ) =U' 2 {a i ,a i \o' i ,o' j ), (149) 

the operators V- and V-j are Hermitian ones, 

N 

H'eff E v'w ( 15 °) 

i=l N>i>j>0 

and operators V- and V-j fulfill one more condition being as follows. Let {y?nl n — 0, 1, ...}be 
an orthonormal basis of eigenfunctions of operator U[ in L 2 (X, dfx) defined by the following 
equation: 

= J U[(a,a'Ma')df,(a'). (151) 

Let U 2 is the operator in L 2 (X x X, d/i x d/i) specified by the relation 

(U' 2 (p)(<Ti, cr 2 ) = J U' 2 (a u a 2 \a[, a^fioi, cr' 2 )dii(a)dii(a'). (152) 

Let {U' 2)mn,m'ri be a matrix element of operator U' 2 with respect to the basis {(Pi®(Pj\i,j = 
0,1,...}. 

{U' 2 )mn,m'n' = (<Pm ® Vn\U' 2 \Vm> ® ¥n>) ■ (153) 
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Then, for any n = 1,2, 3... 

Moo.on = (154) 

and similar equalities take place where index n stands on remaining three places and other 
have zeros. 

Proof. Let Ui be an operator in specified by the following relation: 

(U[<p)(a) = J U[(a,a'Ma')dfi(a'), (155) 

and let {y2„|n = 0, 1, ...} be an orthonormal basis of eigenfunctions of operator Ui. 
We represent the operator as a sum of three summands: 

% = Vff + Kf + (156) 

where operators Vff, V?f, Vff are based on functions (U 2 )(<Ji, a 2 ), (U 2 )(ai, a 2 ), (U 2 )(ai, a 2 ) 
defined below by the principle which the is based on U 2 (ai,a 2 \a[,a' 2 ). The functions 
U 2: U 2 , U 2 are defined by the following condition. Let U 2 , U 2 , U 2 be operators based 
on functions U 2 , U 2 , U 2 using the same scheme as used by the operator U 2 based on the 
function U 2 . Among all matrix elements of operators U 2 , U 2 only the following elements 
are non-zero: (U 2 ) 00 ,om, (U 2 l ) 00tm0 , (C/|) m0)0 o, (U 2 ) 0m ,oo, rn = 1,2,... and among all matrix 
elements of the operator C/f only remained elements are non-zero. It is easily seen that if 
we redefine t^ 6 ' 3 in a proper manner, then we can include summands corresponding to V^j 1 , 
Vff into % 

But Vij and, consequently, are first-order values by the constant of protein interaction 
A. Therefore, in just described substitution, the basis of eigenfunctions of operator U\ 
{(f n \n = 0,1,2,...} passes to the basis of eigenfunctions of redefined operator U[ {ip n \n = 
0, 1, 2, ...} such that <p n — ip n are the first-order of vanishing values by a coupling constant A 
Vn = 0,1,2,.... 

The matrix elements of just redefined V? of the form (U 2 )oo,om m — 1? 2, ... (U 2 is coupled 
with Vfj by the relation we've mentioned above) relating to the basis {<p n } are equal to zero. 
But since are of first-order of vanishing by the interaction constant, then matrix elements 
of the form ([/ 2 )oo,om of redefined (U 2 ) related to the basis {i> n } are of the second order of 
vanishing by the coupling constant. 
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The procedure described above can be continued endlessly by sequentially lowering the 
order of vanishing of those terms we want to omit in the Hamiltonian. 

The proposition may be considered as a proven one even without discussing a problem of 
convergence of the above described iteration procedure, as all formulas which may be derived 
hereafter will be just asymptotic expansions in a small parameter A. 

So, the effective Hamiltonian for considered system of protein molecules, which we want 
to use for defining the ground state of our system, has the following form: 

N 

i=l l<i<j<N 

Here operators Vi have "act" only on coordinate <Ji of the wave function and differ from 
each other only by a number of the coordinate on which they "act". Operators Vij "acts" 
only on coordinates of numbers % and j of the wave function and differ from each other only 
by numbers of coordinates on which they "acts". 

Operator V# "acts" only on coordinates of numbers i and j of the wave function and 
differ from each other only by numbers of coordinates on which they " acts" . 

If in (I157p we evidently extract the dependence from the interaction constant and number 
of particles, then we shall find 

N X 

i=l ±<i<j<N 

where prime of H e ff is omitted and does not depend on iV (asymptotically at N — > oo). 

We intend to investigate the Hamiltonian H e /f using the apparatus of secondary quan- 
tization. 

Let {ip n \n = 0, 1, 2...} be still an orthonormal basis of eigenf unctions of operator U\. Let 
J 7 = r(Z/2pr, d/u,)) be a boson Fock space over Lz(X, dfi). Let af , a« be operators of particles 
creation and annihilation in the state acting in J 7 . Vi = 0, 1, ... operators af and cij are 
conjugated to each other and fulfill the following canonical commutating relations: 

[a h aj] = [af,af] = 0, i,j = 0, 1,2..., 

[a i( a+] = Sij, (159) 

where 5ij is a common Kronecker delta. We suppose that eigenfunctions (p n are numerated in 
such a way that E n increases with increasing of number n. In representation of the secondary 

39 



quantization the Hamiltonian H e ff is given by: 

oo . 

Heff — 2_j E n O'ti a n + T^y a>m a n(U2)mn,m>n'a>m>a n >- (160) 

n=0 m,n,m' ,n' 

But it is a common Hamiltonian studied in the superfluid theory, and we can use the same 
methods for its investigation. 

At T = a portion of proteins is in the ground state ipo. Let N be a number of protein 
molecules in the ground state. As ^ < 1, then N N ^° <C 1. Because of this reason in the 
right part of relation a^aQ — a^ao = 1 we can neglect a unity and suppose and are 
usual c-numbers. 

So, we can simply suppose that ao = \fNe~ %v and = y/Ne tip and for some real number 

cp. 

Now let's consider the second summand in the right part of (1160!) . According to the 
assumption proven above we can suppose that this summand has no terms linear by a$ , a», 
% = 1,2,.... Quadratic terms have an order of A, cubic -^=, and quartic j^. Therefore, 
supposing the number of proteins in the cell is finite but a sufficiently large and interaction 
constant A is small, in the Hamiltonian (I160p we can keep only quadratic terms on operators 
, (X{^ % 1,2, — 

Consequently, the effective Hamiltonian H e ff is given by: 

oo 

H e ff = 71CN + ^ E nan a n 
n=l 

oo oo oo 

+A{ 2j v4 mn e" 2 ^a+a+ + ^2 A * mn e 2l,p a m a n + 2J B mn a^a n }. (161) 

m,n=l m,n=l m,n=l 

Here c is a material constant (does not depend on N ), A mn , B mn are some matrixes A mn = 
A nm , B* mn = B nm , * is a sign of complex conjugation. 

There is a classical variable <p which is canonically conjugated to variable J = 7rN . 
Variable p satisfies the equation <p = c and, therefore, value ip = tp — ct (t is time) is a 
motion integral. Whole our system (denoted by M) has a structure of a direct product of 
two systems M = Mi x M 2 where, roughly speaking, the system Mi is described by variables 
ip and J, and the M 2 system is described by (noncommutative) variables 

5 ^ lj 2 

Now, let's consider two cases: when the cell is a live and when it is dead. In the case of 
the living cell, the expression under the integral for statistical weight contains a multiplier 
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5((f — if'), if' is a fixed real number, and now the free energy F(ip') is a function of if'. 
However, one can readily see that the free energy F(ip') does not depend on if' . Indeed, the 
change of if' to the value Sip can really be compensated by proper canonical transformation 
of operators a+, a n , n = 1,2...: 

a n ^ a n e~ iSv . (162) 

Now, integral if fulfills the equivalence principle we have considered above (Prokhorenko 
and Matveev, 2011). As the formula of the generalized microcanonical distribution includes 
a multiplier factor 5(f — if'), the Mi system may be described by using the Hamiltonian 
(11611) where if is given as a some value and N is given taking into account that the total 
system particles number is equal to N. 

Now, let's consider a dead cell. In this case the cell is described by means of the equilib- 
rium Gibbs distribution. Let's write out Hamilton's equations for ip and J. So, we have: 

f = c, 

j = + \L(a + ,a), (163) 

where L(a + ,a) is a quadratic function of operators a+, a n , n = 1,2.... But if TV — ?► 00, 
then J ~ N, and L(a + ,a) ~ 1. Therefore, when N values are high enough, L(a + ,a) value 
may be neglected in the Hamiltonian equations for if and J and the following equations are 
obtained: 

if = c, 

j = 0. (164) 

Thus, when iV values are high enough, dynamics of the Mi system separates from the 
dynamics of M^. This means that dynamical variables for M\ and M2 systems are indepen- 
dent and the probability distribution for the Mi system, corresponding to Gibbs distribution 
for the whole system, is defined by the formula: 

p(f, J) = const5(J - J'). (165) 

Now, let's define the distribution function for the system M^. For this we use a classical 
analogy. Suppose the Hamiltonian system K has a structure of direct product of systems 

41 



Ki and K 2 , K = K 1 x K 2 and H (x, y) is the Hamiltonian of the whole system, x e K ± , y e 
i^ 2 - If the K system is described by the Gibbs canonical distribution corresponding to T 
temperature, then the probability distribution for the K 2 system is the following: 

1 _ F(y\T) 

p{y) = t , (166) 

where 

F(y|T) = -Tin y dr s e-^, (167) 

and Z is a normalized factor. Alternately, as can be seen above, the effective Hamiltonian 
for the K 2 system can be obtained in the following way. Let's write out the Hamiltonian 
equations for the system K 2 . Assume p and q are canonical coordinates and momenta for 
the system K 2 system. Then: 

. = dH(x,y) 
dq 

i=^. (168) 

Now, if we average right parts of the two last equations by the conditional Gibbs distri- 
bution for K\ w(x\y) having a specified value y G K 2 , then we obtain closed with respect to 
y Hamiltonian equations where F(y\T) is the Hamiltonian. 

Therefore, in order to define the effective dynamics for the M 2 system for a dead cell, 
the following method, reasoning by analogy, can be used. Heisenberg equations for a+, a n , 
n = 1,2... can be written out and averaged over ip and J. However, as shown above, 
if N — > oo, then if is independent of a+, a n , n = 1,2... and distribution for if is uniform. 
Direct calculation shows that if the right parts of Heisenberg equations for a+, a n , n — 1, 2... 
are averaged by the uniformly distributed ip, we obtain Heisenberg equations where the 
Hamiltonian is: 

oo 

H eff = C N + ^2 E ^ a n a n 
n=l 

oo 

+A ^2 B mn a+a n . (169) 

m,n=l 

It's evident, when T = N = Nq, i.e. all protein molecules are folded. 
Now, let's pass to analysis of H e ff using a well-known theory of normal forms of quadratic 
Hamiltonians. For Hamiltonians of a general type this theory was developed by H. Poincare 
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and G.D. Birkhoff (see Arnold, 2003) and for quantum mechanics it was adapted by N.N. 
Bogoliubov (Bogoliubov and Bogoliubov (jr), 1984). This theory was established only for 
dynamical systems with a finite number of degrees of freedom, but we use it for systems with 
an infinite number of degrees of freedom, as it is a standard practice in physics. For this 
purpose we replace the Hamiltonian i? e // by cuted Hamiltonian describing the system 
having L degrees of freedom: 

L 

H c eff = ncN + E E na+a n 

71=1 

L L L 

+A{ A mn e- 2 ^aia+ + E A* mn e 2 ^a m a n + E 5 mn a+a„}. (170) 

m,n=l m,n=l m,n=l 

In this case, the theory of quadratic Hamiltonians confirms that for general matrixes A mn , 
B mn there is such a linear transformation of creation and annihilation operators a+, a n , 
n = 1,2... to operators £+, £ n , 

L L 

a r . 



m=l m=l 
L L 

Zn = J2 U nm a m + E V ^ 
m=l m=l 



'in) 



that Wn = 1,2... £ n is conjugated to this transformation is canonical: 

[£nj Cm] [Cn ' Cm] 0) 

= <U, m,n= 1,2,3... (172) 
and using new variables, Hamiltonian H c e ^ is: 

^// = E W «£ ^ + E + CnCn). (173) 

i=l i=K+l 

u n , Xn are some real numbers. However, the physical considerations make clear that H^sf 
should be low bounded. Let's denote "H c as Hilbertian space where Hamiltonian H^f acts. 
This Hilbertian space is represented by the tensor product of Hilbertian spaces Hf 

L 

^ C = (g)^, (174) 
i=i 
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where every "H? is isomorphic to L 2 (R) and this isomorphism can be chosen in such a way 
that the following condition is fulfilled. In L 2 (R) the operators of coordinate q and momenta 
p operate in the standard way: 

pq — qp = -. (175) 

i 

an d £n operate only on the n-th multiplier in ( 11741) . The isomorphism between and 
L-2,(R) mentioned above can be chosen in such a way that due to this is isomorphism the 
following is correct: 

£t = -^(p + ig). (176) 

But due to the isomorphism between T-L c n and L n (R) mentioned above the operator £%£n + 
£n£n equals to the operator p 2 — q 2 . Evidently, the last operator is not neither low nor upper 
bounded. Therefore, the physical requirement of positiveness of -f^T// results in the fact that 
in the right part of (11731) only the first summand in the right part is non-zero. As a result: 

L 

H c eff = J2^dn, (177) 

i=l 

and for any n — 1,2, 3... u n is a real positive number. 

Operators £+, are creation and annihilation operators of some quasi-particles. As the 
formula (j!77p indicates, at zero temperature all occupation numbers of these quasi-particles 
are equal to zero: n\ — 0, 1 = 1, 2.... 

Now, we are interested in the average: 

(a+a n ), n = l,2... (178) 

for the ground state of our effective Hamiltonian. Evidently, the condition = 0, 1 = 1, 2... 
implies Vra = 1,2... 

oo 

(a+a„) = \V' nm \\ (179) 

m=l 

where matrix V^ m is defined according to following equation: 

L L 

< = Yl U '™Z™ + E y'nmU- (180) 
m=l m=l 
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L 

Clearly for the interaction of general form U 2 between protein molecules Yl \V'nm\ 2 > 0. 

m,n=l 

This means that (a+a n ) > at least for one n = 1, 2.... 

So, if the cell is live, then the number of protein molecules in the unfolded state is 
non-zero. 

As a conclusion of this section we note that according to above reasoning, the number of 
protein molecules in the unfolded state approaches the finite value, if the cell volume tends 
to infinity. The last should mean that the number of the unfolded protein molecules in the 
real cell is negligible. The possible solution of this difficulty is the fact that the protoplasm 
has a quasi-crystalline structure only locally within the range of a small volume we call a 
domain consisting of several physioatoms. Therefore, we can supply the stated analysis only 
within the range of one domain. Since the domain number in the cell is proportional to its 
volume, then the number of protein molecules in the unfolded state is also proportional to 
its volume and has anon-zero value (per unit volume) in the thermodynamic limit. 

7 Discussion of Relation between Mean Field Model 
and Superfluid Bose Gas Model on Protein Config- 
uration Space 

In the previous section we quoted the nonlinear Schrodinger equation describing our system 
in a mean field approximation: 

SeL 2 (X,dfi), (181) 

where £ G L 2 (X,d/j,), and V 1 , V 2 (£) are operators in L 2 (X,d/j,). 

The question arises, to which degree the results obtained within the framework of the 
mean field method (superfluid model) correlate with results obtained within the model of 
the superfluid Bose gas on the protein configuration space. Now, we'll show that these two 
models give the same result for spectrum of single-particle excitations. 

Suppose U 2 used to construct V 2 (£) fulfills the condition which defines, as was said, that 
all the matrix elements of U 2 such that for them three of four indexes equal to zero and 
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fourth is non-zero equal to zero. Here we discussed matrix elements regarding the basis of 
eigenf unctions of operator V 1 . 

Since the interaction constant A is small and almost all protein molecules are in the 
normal unfolded state, then £ is given by: 

i = ifo + V, (182) 

where i] is orthogonal to (fio, </2 , tpi, </?2--- are an orthonormal basis of eigenf unctions of oper- 
ator V 1 , and ifo corresponds to the minimum eigenvalue of V 1 . The value 77 is a first-order 
value by the coupling constant A. Let r] n be Fourier coefficients of vector 77 with regard to 
basis ip n 

00 

V = y^?7n^n- (183) 

n=l 

If in (11821) we keep only terms linear on rj, then complex conjugate to (11821) has the 
following form: 

Q OO OO 

- i— 77* = E nV * + B mn V*m + 2 A ™nV m , Tl = 1, 2, ... (184) 

m=l m=l 

A similar derivation can be made for the equation for 4;T] n . On the other hand in superfluid 
model, after suitable canonical transformation, a+,a„, satisfy the Heisenberg equation: 

a£ = i[H e ff,a£], a n = i[H eff ,a n }. (185) 

If we write down these equations in close detail, then we receive: 

q OO OO 

- = #n<4 + Yl B ^ a m + 2 Y A *mn a m, U = 1, 2, ... (186) 

m=l m=l 

and an equation obtained from the previous one by the Hermitean conjugation. 

But equations (118411 and (1186D have the same form, therefore, if variables {a+,a n } and 
{?7+,?7 n } are subjected to linear transformation of the same form, then equations (I184j) and 
(11861) for transformed variables coincide. Particularly, if the linear transformation transfers 
{a+,a n } into {Cj~,£n}, then r\ n transformed under the same transformation are changed 
according to the following law: 

r] n = conste"^"'. (187) 
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So 7] is given by 

oo 

V = l£fne- iUnt . (188) 

n=l 

This formula shows that Bogolyubov's frequencies also describe a spectrum of single-particle 
excitations in the mean field model. 

8 Comparison of Van der Waals Protoplasm Model 
and Superfluid Bose Gas Model on Configuration 
Space of Protein Molecule 

The present work considers two models of Ling's cell protoplasm microstructure. 

In the first model called Van der Waals model, we supposed that nontrivial first integrals 
of the system are so that their fixing by predefined values characterizing the resting state 
of a cell leads to the fact that protein molecules are situated in points of a lattice in the 
unfolded conformation. In addition, if the rapid descending of Van der Waals interaction 
between protein molecules with distance increasing is taken into account, the thermodynamic 
features of the living protoplasm can be calculated just as for the ideal gas of proteins. The 
thermodynamic features of the dead protoplasm can be calculated using a well-known Van 
der Waals interpolation formula for the free energy of a system. With this model we've 
obtained an expression for the quantity of heat the cell released while dying, and for a 
number of potassium ions releasing from the cell during this process. 

In the second model we emphasize the analysis of internal protein molecules structure. 
As we have already mentioned, the protein molecules are supposed to be situated in points 
of a lattice, but have nontrivial internal degrees of freedom; and we study the structure of 
the ground state of this molecular system. The present work shows that within the range 
of weak interaction between protein molecules, the ground state of the interacting proteins 
system should be outlined as a ground state of the Hamiltonian of the Bose gas with a weak 
interaction on the protein molecule configuration space. To analyze this Hamiltonian we 
used standard methods of the superfluid theory. Taking into account that the interaction 
between protein molecules is weak and the bulk of protein molecules of the system is in the 
ground folded state, our Hamiltonian can be replaced with an effective quadratic one, but 
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it has different forms depending on if the cell is live or dead. In addition, it clears up which 
nontrivial motion integrals in the involution should be included to the generalized Gibbs 
distribution describing the living cell. The use of these effective quadratic Hamiltonians 
reveals that in the dead cell all the protein molecules significant for model properties are 
folded, while in the resting living cell the number of unfolded protein molecules is non-zero. 

The question arises: what is the relation between considered protoplasm models? Partic- 
ularly, how do conservation laws postulated for them link? We think, the considered models 
complement each other and connection of the named conservation laws comes down to the 
following. If we fix the motion integrals for the second model by some predefined values, 
then protein molecules are in the unfolded state, and we can suppose that the nature of 
interaction between different protein molecules provides the formation of a lattice being an 
energetically favorable structure. In other words, the conservation laws for the first model 
appear to be consequences of the conservation laws for the second model. If we conversely 
consider the Van der Waals protoplasm model as the main one and assign conservation laws 
to it, according to which the protein molecules in the living state are situated in points of a 
lattice, then it is reasonable to suggest that such a space configuration of protein molecules 
would be favorable for the unfolded (not folded) state of protein molecules for a number 
of reasons. In other words, the conservation laws for the second model are consequences 
of the conservation laws for the first model. This situation is certainly very inexact and 
hypothetical. 

The questions, on which principles should new models be constructed (except ones con- 
sidered here) and how to develop the kinetic theory of the living cell, are the subject of 
further investigations. 

9 Conclusion 

In the present work we have constructed and investigated properties of two complementary 
models of protoplasm — physical basis of life. The work was realized basing on generalized 
thermodynamics we proposed (Prokhorenko and Matveev, 2011). Within the framework of 
stated assumptions we explained a number of properties and phenomena observed in living 
cells, formalized (in the extended sense) by the physical theory of the living cell by Ling 
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(2001). However, the results we had formerly obtained, were inequalities and denoted only 
processes direction (in excitation and death of a cell heat releases, volume reduces and so on), 
there were no certain numerical evaluations obtained. To obtain them we need to specify 
the properties of intermolecular interactions which escape from the point of view of available 
analyze methods. That's why, we should construct different protoplasm models emphasizing 
one system parameter after another. 

In the present work we constructed and investigated two models of protoplasm: one 
of them is naturally called Van der Waals model, and other is the model of superfluid 
Bose gas on the protein molecule configuration space. Our aim in this work was not to 
construct a model giving the most exact agreement with experimental data but to show that 
the constructing of such models is reasonable and possible. Qualitative agreement of the 
obtained results with experimental data gives an evidence of vitality of the thermodynamic 
theory of the living cell we proposed (Prokhorenko and Matveev, 2011). 

Our theory can face the misapprehension from devotees of the non-equilibrium thermo- 
dynamics, as it is based on the equilibrium statistical physics and thermodynamics. Their 
main argument is evident: in the equilibrium state the maximum entropy is reached, there- 
fore, the cell cannot perform the biological work. However, the other fact is also evident: 
the non-equilibrium thermodynamics still haven't given the quantitative description even for 
elementary phenomena, for example, electric resting potential. In addition, if we accept that 
the cell is living under the laws of non-equilibrium thermodynamics, we must also recognize 
that 200 years experience in successful modeling properties of the cell by non-living systems 
is untenable in principle. It would ruin all our knowledge about the living cell without giv- 
ing anything in return. As for our generalized thermodynamics, the essential clarification is 
required: though the states which it operates are static by time, they are not equilibrium in 
the sense that their entropy is not the maximum one of all the states having the same energy 
(this is an essential feature of the state of the living material). The existence of such states 
even for the most realistic statistical mechanics systems is proven by one of us (Prokhorenko, 
2009). 

We are very grateful to P. Agutter, A.V. Koshelkin, Yu. E. Lozovik, A.V. Zayakin, E.N. 
Telshevskiy and N.V. Puzyrnikova for valuable critical comments on this article and very 
useful discussions. 
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10 Appendix 1. Transition of Cell from Living to Dead 
State in Weak External Alternating Field 

Here and in our previous work (Prokhorenko and Matveev, 2011) we considered two extreme 
states of a cell: non-equilibrium state of rest and equilibrium state corresponding to death. 
Now, let's ask a question, can we describe death of the protoplasm within the framework 
of our generalized thermodynamics, i.e. to construct an example of a process transferring 
a cell from the living state of rest to the state corresponding to death. If we could not 
construct an example of such a process, all our theory would appear to be doubtful. In 
addition, up till now we considered only equilibrium statistical mechanics of the living cell, 
and the constructing of an example of the transformation process is the first step toward the 
construction of the non-equilibrium statistical mechanics of the living cell. 

Let's start the analysis of the cell death mechanism from a certain problem of the non- 
equilibrium statistical mechanics. Let's consider the Hamilton system having n degrees 
of freedom which is described by the Hamiltonian H(p, q) in the field of external forces 
ef(t), eeK such that the complete Hamiltonian of the system is as follows: 

r = H + ef(t)P, (189) 

where P(p, q) is a function of canonically conjugated momenta and coordinates of a system. 
We suppose that applied force fit) is a real function of time of the following form: 

f(t) = ^2a v coa(yt + <p v ), (190) 

where phases y?„ are independent random values uniformly distributed by a circle. We 
suppose the frequency spectrum of the applied force is essentially continuous, so that sums 
of the following form 

5>(i/)a* (191) 
with continuous F(y) in the limit can be replaced by integrals 

+oo 

J F{v)l{v)dv. (192) 



This problem was considered by N.N. Bogoliubov (1945) and here we briefly quote the results 
obtained by them (relevant only to classical mechanics). 
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Let's denote D t is a probability density of coordinates and momenta in a time t provided 
that phases (p v have some defined values. The probability density for distribution p and q 
in the common sense, i.e. when phase values ip u are inessential, can be found from D t by 
averaging over all phases: 

Pt = A- (193) 

In the initial moment t = we suppose that distribution of coordinates and momenta does 
not depend on phases, in other terms 

D = po- (194) 

Temporal evolution of D t should be passed in accordance with a well-known Liouville's 
equation: 

dD 

-^l = (T,D t )+ef(t)(P,D t ), (195) 
where (A, B) is a Poisson bracket defined by the following formula: 

Let's introduce more one-parameter group of operators T t t e K acting on dynamical 
variables according to the formula: 

T t F(p,q) = F(p t ,q t ), (197) 

where p t , q t is a solution of canonical Hamilton equations 

dpt _ dH(p, q) 
dt dq 

dq dH(p, q) 

dt dq 

having initial data 



(198) 



Po =P, % = Q- (199) 

Using these assumptions and designations (Bogoliubov, 1945), in the limit of small e the 
following equation for p t was derived: 

t 



dpt 
dt 



= [H, Pt ]+e 2 J A(t-r)(P,(T T _ t (P,p T ))dr, (200) 



o 
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where 

+00 

A(r) = ^ J I{u) cos(z/r)rfr. (201) 







Further Bogoliubov (1945) considered the case when the system exposed to external force is 
a harmonic oscillator of n-dimensions having incommensurate frequencies uj = (oji, u> n ). 
Let's denote through a o t = ~pi function of action variables obtained from p t by averaging 
by angular variables 9±, 6 n . Function P(p, q) becomes a function of action-angle variables 
and can be expanded to the Fourier series: 

P = Pne^ H . (202) 

a 

In this case, at the same initial conditions on D t , in the limit of small e the evolution of a t 
is described by the Fokker-Planck equation: 

k=l s=l s K 

where 

A sk (I) = ^n fc n s /(n-a;)|P^(/)| 2 . (204) 



This equation, particularly, describes diffusion of distributions which at the initial moment 
can be distributions concentrated in one point, i.e. 

n 

^ > fc=l 

But the last distribution is a generalized microcanonical distribution which we used to de- 
scribe the states of rest of the Ling's cell. Therefore, it seems to be natural if the method 
proposed by Bogoliubov (1945) can be used to construct an example of a cell transformation 
process from the living to dead (equilibrium) state. Now, we shall show how it can be done 
using the results of Appendix 1 of the work by (Prokhorenko and Matveev, 2011). 

So, assume there is a Hamilton system M with n + k degrees of freedom where k <^ n, 
which is described by the Hamiltonian H and where k independent first integrals in the 
involution are defined. There was shown in the Appendix 1 of the above-mentioned work 
that some covering Hamiltonian system M' of system M can be represented as a direct 
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product of M' = M[ x M 2 , so that the canonically conjugated coordinates on M 2 are first 
integrals K[, K' k (more properly, their lifting to M') playing a role of momenta and "angle" 
variables ipi,..., ipk, passing all the real axis and playing a role of coordinates. 

Let 7r be a canonical projection of M' on M. The function / on M' we (Prokhorenko 
and Matveev, 2011) called a periodic one if f(x) = h(%(x)) for some function h defined 
on M. It has been suggested (Prokhorenko and Matveev, 2011) to represent the mixed 
state of the M' system using some positive periodic function p on M'. Put by definition 
Rl '■= {% £ ^'Hvil < L, \ifik\ < L}. If A is a periodic function on M', then its average 
over the state corresponding to distribution p was suggested (Prokhorenko and Matveev, 
2011) to calculate using the following formula: 

/ p(Ki, K k , (p t , (p k )A(Ki, K k , (pt, (p k )dK 1 ...d(p k 
l-»oo J p(K 1 ,...,K k ,(p 1 ,...,(p k )dK 1 ...d(p k 

Rl 

In this work we have shown that the limit (12061) always exists. We shall call the distribution 
function p(K 1 , K k , <pi, fk) as a normalized one if 

lim / p(K u ...,K k ,<p-i_, ...,<p k )dK 1 ...d<p k = 1. (207) 

Rl 

The entropy of a state corresponding to the normalized distribution p{K\, K k , (p 1: ip k ) 
is defined by the following formula 

S = - lim — / p(K 1 ,...,K k ,ipi,...,(p k )lnp(K 1 ,...,K k ,ipi,...,ip k )dK 1 ...d(p k . (208) 

Rl^ck L k J 
Rl 

The last limit exists while the proof of its existence is the same as for the limit 



The Hamiltonian of the covering system M' we also (without a risk to make an error) 
denote as H. If our Hamilton system is a Ling's cell, then n ^> k. In this CclSC, clS shown 
(Prokhorenko and Matveev, 2011), the dynamics of the system M' 2 can be considered sep- 
arately. This dynamics is a Hamiltonian one and is defined by Hamiltonian F(x\T) where 
F(x\T) is the free energy of the system M[ at temperature T provided that the M' 2 system 
is situated in point x: 

F(x\T) = -Tin J rfrje-^, (209) 

y G M[, x G M 2 , dTy is an element of phase volume on M[. 
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But Ki,...,Kk are the motion integrals, therefore F(x\T), does not depend on angular 
variables <pi, but depends only on K\, K k . Put by definition 

F'{K x ,...,K k \T):=F{x\T). (210) 

Further, we are interested in the case when F'(Ki, K k ) reaches its minimum in the whole 
area (of non-zero volume) O. In the work by Prokhorenko and Matveev (2011) the cor- 
responding thesis is called the equivalence principle and there was shown that only under 
condition of F'(K\, ...,Kk) constancy in the whole area O of the non-zero volume our gen- 
eralized thermodynamics gives new results in comparison with the common one and can 
provide the thermodynamic descriptions of resting state of the Ling's cell. 

So, F(x\T) is a constant on the direct product 0xl fe and, therefore, it defines a trivial 
dynamics in this area, i.e. if (p , q ) G O x R fc , then Vt G R (p t , It) = (po, Qo)- 

Using the last circumstance the equation (1200 j) can be simplified as follows 

t 

9 -§ = e 2 f A(t-r)(P,(P,p T ))dr. (211) 
o 

However, in the limit of small e p T {z) varies with time slowly and is almost constant wherever 
A(t — r) noticeably differs from zero. Therefore, in the integral from the right part of the 
last equation we can replace p T by its value at r = t. As 

+oo 

/ A(t)dt = *-f> (212, 



we can find 

to = *m {P , {P , pd) . (2i3) 

This equation of Fokker-Planck type is correct in the area O x M. k . 

As for behavior of function p t outside the area Oxtf, the physical considerations make 
clear that the probability to find the system M' outside the domain O x IR fc is negligible. 
Therefore, it is sufficient to investigate the behavior of p t inside the domain O x M fc by 
applying the appropriate boundary conditions for p t on the boundary of O x R fc which 
provide the self-adjointness of equation (12131) and keeping the probability: 

J Pt (x)dTl = l, (214) 

OxR k 
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where dT 2 x is an element of phase volume on M' 2 . However, we won't make an ascertaining 
of the form of these boundary conditions but just make an assumption that the domain O 
is a compact manifold (without boundary). 

The self-adjointness condition of operator from the right part of the Fokker-Plank equa- 
tion is a consequence of equality of forward and backward probabilities of arbitrary transition. 
This equality of forward and backward probabilities of transitions is a consequence of com- 
plete Hamiltonian invariance with respect to the time sign conversion operation (See the 
principle of kinetic coefficients symmetry by L. Onsager). 

Now, following Bogolyubov (1945), we can show that the entropy of a state defined 
by relation ( 1208 j) is increased monotonically in course of time due to equation ( I213p . The 
entropy ( I2U8P reaches its maximum at the constant distribution function, i.e. at t = +00 

p+aoiKx, <p k ) = const. (215) 

But if the distribution function p +OQ (K 1 , ip k ) is constant, then it corresponds to the 
Gibbs equilibrium microcanonical distribution. 

So, an example of the Ling's cell transformation from the living (resting) to the dead 
state is its evolution in the weak external field which is a sum of harmonic oscillations with 
continuous frequency spectrum and random independent phases uniformly distributed along 
the circle; at that, for spectral density of intensity of this field I(v), 1(0) 7^ is required to 
be fulfilled. 

The fact that the Fokker-Planck equation includes the spectral density of intensity only 
through 1(0) is very significant and can be verified experimentally. This conclusion is in a 
good agreement with observed destructive influence of infrasound, magnetic storms and any 
white noises to living organisms. 

11 Appendix 2. On ATP Structuring Role as Part of 
(ATP)m-(Protein)n-(H20)p-(K+)q Complex 

In this Appendix we discuss the theoretical relation between the capability of an ATP 
molecule to define the structure of a physioatom in the resting state and non-ergodicity 
of the Ling's cell. The purpose of this Appendix is to represent a demonstrative physical 
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view of the structuring ATP influence on the physioatom. We do not give a mathematically 
consistent characteristic of ATP-protein-water-K + interaction here. 

We are interested in the following issues. How much is a number of protein molecules in 
the water-protein complex managed by one ATP molecule? Why the physical disturbance 
generated by the ATP propagates without dissipation to the large number of proteins and 
how does this disturbance influence on their conformation? 

The fact that, according to our main assumption, the Ling's cell represents a Hamiltonian 
system having a vast number of first integrals in the involution makes us think, can we use 
the theory of completely integrable (in the sense of Liouville) systems to describe such a cell? 
The Korteweg — de Vries equation (Arnold, 2003) which can be an example of a completely 
integrable system 

u t = 6uu x — u xxx (216) 

was originated in the shallow water theory (in narrow ship channels). This equation is re- 
markable by allowing the solutions in a form of solitary waves (solitons) propagating without 
dissipation. The only parameter characterizing the soliton is its velocity. The Korteweg — de 
Vries equation also allows for multi-solitonic solutions which break into separate solitons 
propagating with different velocities at t — > ±00. A significant property of multi-solitonic 
solutions is the fact that in case of solitons collision their velocities do not change. 

We consider that the distribution of the ATP molecule physical influence on surrounding 
protein molecules has something similar to the propagating of solitons because there are 
many commutative first integrals for the Ling's cell, therefore, this case should be somewhat 
similar to the case arising in the theory of completely integrable systems. 

So that one ATP molecule could effectively manage the surrounding complex of water and 
proteins, the disturbance transferred by solitons should not dissipate. For analysis simplicity 
let's suppose that the distribution of the physical impulse from the ATP is described by a 
one-dimensional Korteweg — de Vries equation. Then every soliton is unambiguously defined 
by the only parameter, its velocity. Therefore, to prevent loss of information from the ATP 
molecule, solitons' velocities shall not change after their collision. But the last property is 
fulfilled for multi-solitonic solutions of the Korteweg — de Vries equation, as we've mentioned 
above. 

The theory of completely integrable systems is a rather developed one (Bullaf and Caudry, 
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1999). A lot of integrable equations were constructed on the line, for example: nonlinear 
Schrodinger equation, sine-Gordon equation, Toda chain and so on. The main method of 
integrating such equations is a method of inverse scattering problem (Bullaf and Caudry, 
1999). The property of isolated solutions to pass through each other without velocity changes 
can be common for all of them and is explained by the presence of the complete set of 
independent commutative first integrals. 

Let's take, for instance, the finite Toda chain consisting of N particles in the line (Mozer, 
1975). The state of this system is fully described by defining N particle coordinates {x^i = 
1, N} and iV momenta {pi\i = 1, N}. By definition the Hamiltonian of this system is: 

N 2 N ~ 1 

H = ^ Pi. + ^ e Xi+1 - Xi . (217) 

i=l i=l 

As shown by Mozer (1975), at t — > ±oo the distances between different particles tend 
to infinity. This system is completely integrable (Mozer, 1975) and there is a set of N 
independent first integrals in the involution for this system. If distances between particles 
are so large so their interaction can be neglected, then those integrals are just elementary 
symmetrical polynomials of momenta (velocities) (Mozer, 1975). 

h=Pi+P2 + — +Pn, 
h = P1P2 + ••• + PiPn + •• + Pn-iPn 



h = P1P2-PN (218) 

As before and after particles' collisions the values of integrals Ii,...,I N coincide, to define 
the velocities after collision we have a system of iV algebraic equations which implies that 
velocities of particles after collision are the roots of algebraic equation 

(y-v 1 )...(y-v N )=0, (219) 

where vi, ...,vn are velocities of particles before collision and v is an unknown variable. So, 
velocities of particles after collision coincide with velocities of particles before collision with 
an accuracy of transmutation. 

Yet notice that for integrable systems with a number of freedom degrees N — > 00 the 
equivalence principle is fulfilled. Suppose Ii,...,I N are action variables, and (pi,...,(pN are 
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angle variables conjugated to them. We can take J 2 , ...,In as independent first integrals in 
the involution, which are included in the generalized microcanonical distribution. We have: 

r N 

S(E, I 2 , I N ) = In / dhd Vl J] dtpiSiHih, I N ) - E). (220) 

J i=2 
Integrating this formula by dl\ results in: 

S(E, I 2 , I N ) = In l dH(I l^ lNh +N\n 2tt. (221) 
I dh I 

But In) is a frequency corresponding to ipi] and it is reasonable to restrict our 

choice to systems for which u\ is an asymptotically constant at — > oo. So: 

S(E, I 2 , .., I N ) = - In co^h, I N ) + N In 2vr (222) 

But since uji(Ii, In) is an asymptotically constant, in the limit of iV — > oo can be ne- 
glected, and the equivalence principle is fulfilled for integrals I 2 , In m the limit of — > oo. 
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